A physical template family for gravitational waves from precessing binaries of 
spinning compact objects: Application to single-spin binaries 
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The detection of the gravitational waves (GWs) emitted by precessing binaries of spinning compact 
objects is complicated by the large number of parameters (such as the magnitudes and initial 
directions of the spins, and the position and orientation of the binary with respect to the detector) 
that are required to model accurately the precession-induced modulations of the GW signal. In this 
paper we describe a fast matched-filtering search scheme for precessing binaries, and we adopt the 
physical template family proposed by Buonanno, Chen, and Vallisneri [Phys. Rev. D 67, 104025 
(2003)] for ground-based interferometers. This family provides essentially exact waveforms, written 
directly in terms of the physical parameters, for binaries with a single significant spin, and for 
which the observed GW signal is emitted during the phase of adiabatic inspiral (for LIGO-I and 
VIRGO, this corresponds to a total mass M < 15M©). We show how the detection statistic can 
be maximized automatically over all the parameters (including the position and orientation of the 
binary with respect to the detector), except four (the two masses, the magnitude of the single spin, 
and the opening angle between the spin and the orbital angular momentum), so the template bank 
used in the search is only four-dimensional; this technique is relevant also to the searches for GW 
from extreme-mass-ratio inspirals and supermassive blackhole inspirals to be performed using the 
space-borne detector LISA. Using the LIGO-I design sensitivity, we compute the detection threshold 
(~ 10) required for a false-alarm probability of 10 _3 /year, and the number of templates (~ 76,000) 
required for a minimum match of 0.97, for the mass range (mi, m-i) = [7, 12]M© x [1, 3]Mq. 

PACS numbers: 04.30.Db, x04.25.Nx, 04.80.Nn, 95.55.Ym 



I. INTRODUCTION 



Binaries consisting of a black hole (BH) in combi- 
nation with another BH or with a neutron star (NS) 
are among the most promising gravitational-wave (GW) 
sources for first-generation laser-interferometer GW de- 
tectors such as LIGO 0,0, VIRGO @, GEO 600 @,|| 
and TAMA 300 5]. For LIGO-I and VIRGO, and for bi- 
naries with total mass M < 15M Q , the observed GW 
signal is emitted during the adiabatic-inspiral regime, 
where post-Newtonian (PN) calculations can be used to 
describe the dynamics of the binary and predict the grav- 
itational waveforms emitted [SOUS El- 

Very little is known about the statistical distribution 
of BH spin magnitudes in binaries: the spins could very 
well be large, with a significant impact on both binary 
dynamics and gravitational waveforms. On the contrary, 
it is generally believed that NS spins will be small in the 
NS-BH and NS-NS binaries that are likely to be observed 
with first-generation GW detectors. For example, the 
observed NS-NS binary pulsars have rather small spin, 
Sns / to ns ~ 10~ 3 6]. One reason the NSs in binaries 
of interest for GW detectors should carry small spin is 
that they are old enough to have spun down considerably 
(even if they once had spins comparable to the theoretical 
upper limits, Sns/^ns ~ 0.6-0.7 [lOj, where m^s is the 
NS mass, and where we set G = c = 1), and because 
dynamical evolution cannot spin them up significantly 



(even during the final phase of inspiral when tidal torques 
become important |llj|). 

Population-synthesis studies 01 suggest that in NS- 
BH binaries there is a significant possibility for the BH 
spin to be substantially misaligned with the orbital an- 
gular momentum of the binary. Early investigations 
|l3l Il4| showed that when this is the case and the BH 
spin is large, the evolution of the GW phase and ampli- 
tude during the adiabatic inspiral is significantly affected 
by spin-induced modulations. While reliable templates 
for precessing binaries should include these modulational 
effects, performing GW searches with template families 
that include all the prima facie relevant parameters (the 
masses, the spins, the angles that describe the relative 
orientations of detector and binary, and the direction of 
propagation of GWs to the detector) is extremely com- 
putationally intensive. 

Several authors have explored this issue, and they have 
proposed detection template families (DTFs) that de- 
pend on fewer parameters and that can still reproduce 
well the expected physical signals. An interesting sug- 
gestion, built on the results obtained in Ref. [131 ], came 
from Apostolatos [T^J , who introduced a modulational si- 
nusoidal term (the Apostolatos ansatz) in the frequency- 
domain phase of the templates to capture the effects 
of precession. This suggestion was tested further by 
Grandclement, Kalogera and Vecchio ^S]. The result- 
ing template family has significantly fewer parameters, 
but its computational requirements are still very high, 
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and its signal-fitting performance is not very satisfactory; 
Grandclement and Kalogera 0] subsequently suggested 
a modified family of spiky templates that fit the signals 
better. 

After investigating the dynamics of precessing binaries, 
Buonanno, Chen and Vallisneri henceforth BCV2] 
proposed a new convention for quadrupolar GW emission 
in such binaries, whereby the oscillatory effects of preces- 
sion are isolated in the evolution of the GW polarization 
tensors. As a result, the response of the detector to the 
GWs can be written as the product of a carrier signal and 
a modulational correction, which can be handled using an 
extension of the Apostolatos ansatz. On the basis of these 
observations, BCV2 built a modulated frequency-domain 
DTF that, for maximal spins, yields average fitting fac- 
tors (FF, see Sec. VIB of BCV2) of ~ 0.97 for (7 + 5)M Q 
BH-BH binaries, and ~ 0.93 for (10 + 1.4)M NS-BH 
binaries (see also Tab. VIII, Tab. IX, and Fig. 14 of 
BCV2). Note that the stationary-phase-approximation 
(SPA) templates developed for nonspinning binaries give 
much lower FFs of ~ 0.90 for (7+5)M© BH-BH binaries, 
and ~ 0.78 for (10 + 1.4)M NS-BH binaries, while ac- 
cording to our computations the Apostolatos templates 
give FF ~ 0.81 for (10 + 1.4)M NS-BH binaries. 1 

An important feature of the BCV2 templates is that 
their mathematical structure allows an automatic search 
over several of the modulational parameters (in strict 
analogy to the automatic search over initial orbital phase 
in GW searches for nonspinning binaries), reducing sig- 
nificantly the number of templates in the search banks, 
and therefore the computational cost. However, since 
many more signal shapes are effectively (if implicitly) 
tested against the detector output, the detection thresh- 
old for this DTF should be set higher than those for 
simpler families (for the same false-alarm probability). 
According to simple false-alarm computations performed 
with Gaussian, stationary detector noise (see BCV2) for a 
single template, the gain in FF is larger than the increase 
in the threshold only for binaries (such as NS-BH bina- 
ries) with low symmetric mass ratios m\mij[m,\ + TO2) 2 ; 
while the opposite is true for high mass ratios. [Ulti- 
mately, the issue of FF gain versus threshold increase 
will be settled only after constructing the mismatch met- 
ric for this template family and performing Monte Carlo 
analyses of false-alarm statistics for the entire template 
bank under realistic detector noise.] Although the im- 
provement in FF with the BCV2 DTF is relevant, it is 
still not completely satisfactory, because it translates to 
a loss of ~ 20% in detection rate (for the maximal-spin 
case) with respect to a perfect template bank (the loss 
will be higher if the higher required threshold is taken 



1 The authors of Refs. |15L III! do not include a Thomas preces- 
sion term in the physical model used to test t he temp lates; for 
this reason, the fitting factors quoted in Refs. Il5t are sub- 
stantially lower than our result. T hose authors are currently 
investigating the effect of that term llSl . 



into account). Current estimates of binary- inspiral event 
rates within the distance accessible to first-generation 
GW interferometers hovers around one event per year, 
so a reduction of ~ 20% in the detection rate may not be 
acceptable. 

BCV2 also proposed, but did not test, a new promis- 
ing family of physical templates (i. e., templates that are 
exact within the approximations made to write the PN 
equations) for binaries where only one of the two com- 
pact bodies carries a significant spin. This family has 
two remarkable advantages: (i) it consists only of the 
physical waveforms predicted by the PN equations in the 
adiabatic limit, so it does not raise the detection thresh- 
old unnecessarily by including unphysical templates, as 
the BCV2 DTF did; (ii) all the template parameters ex- 
cept four are extrinsic: that is, they can be searched over 
semi- algebraically without having to compute all of the 
corresponding waveforms. 

In this paper we describe a data-analysis scheme that 
employs this family, and we estimate the number of tem- 
plates required for a NS-BH search with LIGO-I: we as- 
sume 1 M < to ns < 3M Q , and 7M Q < to B h < 12M 
(see Sec. IIII D"|) . In a companion paper [l9|. we show how 
a simple extension of this template family can be used to 
search for the GWs emitted by binaries when both com- 
pact bodies have significant spins (and where of course 
the adiabatic limit of the PN equations is still valid). The 
problem of estimating the parameters of the binaries is 
examined in a forthcoming paper [20(. 

This paper is organized as follows. In Sec. |H]we re- 
view the formalism of matched-filtering GW detection, 
and we establish some notation. In Sec. IIIII we review 
the PN dynamics and GW generation in single-spin bi- 
naries, and we discuss the accuracy of the resulting wave- 
forms, indicating the range of masses to which our phys- 
ical template family can be applied. In Sec. IIVI we de- 
scribe the parametrization of the templates, and we dis- 
cuss the semialgebraic maximization of signal-template 
correlations with respect to the extrinsic parameters. In 
Sec. [V] we describe and test a fast two-stage detection 
scheme that employs the templates, and we discuss its 
false-alarm statistics. In Sec. I VII we build the template 
mismatch metric, and we evaluate the number of tem- 
plates required for an actual GW search. Finally, in Sec. 
IV 111 we summarize our conclusions. 



II. A BRIEF REFRESHER ON 
MATCHED-FILTERING GW DETECTION 

We refer the reader to Ref. H3 (henceforth BCV1), 
for a self-contained discussion of matched-filtering tech- 
niques for GW detection, which includes all relevant bib- 
liographic references. In this section we shall be con- 
tent with introducing cursorily the quantities and sym- 
bols used throughout this paper. 

Matched filtering mElHIHEilillillMIM 

mUlllllslllil is the standard method to detect GW 
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signals of known shape, whereby we compare the detec- 
tor output with templates that approximate closely the 
signals expected from a given class of sources, for a va- 
riety of source parameters. The goodness of fit between 
the template h(X A ) (where X A denotes all the source pa- 
rameters) and the real GW signal s is quantified by the 
overlap 



p[s,h(X A )} 



(s, h(X A )) 



^(h(X A ),h(X A )) 



(1) 



[also known as the signal-to-noise ratio after filtering s 
by h(X A )], where the inner product (g(t), h(t)) of two 
rea l sig nals with Fourier transforms §(/), h(f) is given 
by M 



(9,h) = 2 f 

J — ( 



Sn(\f\) 



*= 4Re X ~Wr dh 

(2) 

throughout this paper we adopt the LIGO-I one-sided 
noise power spectral density S n given by Eq. (28) of 
BCV1. Except where otherwise noted, we shall always 
consider normalized templates h (where the hat denotes 
normalization), for which (h(X A ) 7 h(X A )) = 1, so we can 
drop the denominator of Eq. Q . 

A large overlap between a given stretch of detector 
output and a particular template implies that there is a 
high probability that a GW signal similar to the template 
is actually present in the output, and is not being merely 
simulated by noise alone. Therefore the overlap can be 
used as a detection statistic: we may claim a detection 
if the overlap rises above a detection threshold p* , which 
is set, on the basis of a characterization of the noise, in 
such a way that false alarms are sufficiently unlikely. 

The maximum (optimal) overlap that can be achieved 
for the signal s is y/ (s, s) (the optimal signal-to-noise ra- 
tio) , which is achieved by a perfect (normalized) template 
h = s/ \J (s, s). In practice, however, this value will not 
be reached, for two distinct reasons. First, the template 
family {h(X A )} might not contain a faithful representa- 
tion of the physical signal w. The fraction of the theoret- 
ical maximum overlap that is recovered by the template 
family is quantified by the fitting factor 



FF = 



max\A (w, h(X A )) 



(3) 



Second, in practice we will usually not be able to use a 
continuous template family {h(X A )}, but instead we will 
have to settle with a discretized template bank {h(X A k ^)} : 
where (k) indexes a finite lattice in parameter space; so 
the best template to match a given physical signal will 
have to be replaced by a nearby template in the bank. 
[As we shall see in Sec. II VI there is a partial exception 
to this rule: we can take into account all possible val- 
ues of certain parameters, known as extrinsic parameters 
[22I |26| , without actually laying down templates in the 



bank along that parameter direction.] The fraction of the 
optimal overlap that is recovered by the template bank, 
in the wor st p ossible case, is quantified by the minimum 
match |26ll30| . Assuming that the physical signal belongs 
to the continuous template family {h(X A )}, the minimum 
match is equal to 



MM = minmax(/i(A' ), h(X A k) )). 



(4) 



The required minimum match MM sets the allowable 
coarseness of the template bank [22, H3, : the closer 
to one the MM, the closer to one another the templates 
will need to be laid down. In Sec. I VII we shall use a 
notion of metric [2H I28L IsHI ] in parameter space to char- 
acterize the size and the geometry of the template bank 
corresponding to a given MM. 



III. ADIABATIC POST-NEWTONIAN MODEL 
FOR SINGLE-SPIN BINARY INSPIRALS 

In this section we discuss PN adiabatic dynamics and 
GW generation for NS-BH and BH-BH binaries. Specif- 
ically, in Sees. IIII Al llTTCl we review the PN equations 
and the GW emission formalism developed in BCV2. In 
Sec. ITTTDl we extend that analysis to study the accu- 
racy of the waveforms, and we determine the mass range 
where the waveforms produced by adiabatic models can 
be considered accurate for the purpose of GW detection. 
In Sec. IIII El we investigate the effects of quadrupole- 
monopole interactions (tidal torques) on the waveforms. 
In this paper we restrict our analysis to binaries where 
only one body has significant spin, leaving a similar study 
of generic binaries to a companion paper [2(| ■ As a fur- 
ther restriction, we consider only binaries in circular or- 
bits, assuming that they have already been circularized 
by radiation reaction as they enter the frequency band of 
ground-based GW detectors. 

For all binaries, we denote the total mass by M = m\ + 
7712 and the symmetric mass ratio by 77 = mi 7712 /M 2 ; we 
also assume that the heavier body (with mass m\ > m%) 
carries the spin S\ = Xi TO i; with < xi < 1 (here and 
throughout this paper we set G = c = 1). 



A. The PN dynamical evolution 

In the adiabatic approach |3r3l3^ l?3| to the evolution 
of spinning binaries, one builds a sequence of precessing 
(due to spin effects) and shrinking (due to radiation re- 
action) circular orbits. The orbital frequency increases 
as the orbit shrinks. The timescales of the precession 
and shrinkage are both long compared to the orbital pe- 
riod (this is the adiabatic condition), until the very late 
stage of binary evolution. [These orbits are sometimes 
also called spherical orbits, since they reside on a sphere 
with slowly shrinking radius.] 
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The radiation-reaction-induced evolution of frequency 
can be calculated by using the energy-balance equation, 

where £ is the orbital-energy function, and T the GW 
energy-flux (or luminosity) function. Both have been cal- 



culated as functions of the orbital frequency using PN- 
expansion techniques, and are determined up to 3.5PN 
order 0, 0, 0] ; howeve r, sp in effects have been calculated 
only up to 2PN order [36| • The resulting evolution equa- 
tion for ui, obtained by inserting the PN expansions of 
£ and T into Eq. JSJ and reexpanding [every (Mw) 4//3 
counts as 1 PN] is 



= — r?(Mcj) 5 / 3 
oj^ 5 



1 



743 + 924 77 



336 

34103 13 661 



W 2/3 - ( Y 2 



59 



Xi{t N -S x ) (113^ + 757, 



- 4tt (Mw) 



18144 2 016 

16 447 322 263 
139 708 800 ' 



18 
1712 



V+tzV") (Mw) 4/3 - (4 159 + 15 876 77) tt (Mw) 



1 
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(Mw) 



4 415 358 675 



4032 6 048 



91495 
1512 



■q 2 tt (Muj) 



7/3 



(6) 



where — 0.577 ... is Euler's constant. We denote by 
Ljy oc r x v the unit vector along the orbital angular 
momentum, where r and v are the two-body center-of- 
mass radial separation and relative velocity, respectively. 
Ljv is also the unit normal vector to the orbital plane. 
Throughout this paper we shall always use hats to denote 
unit vectors. (Note for v3 of this paper on gr-qc: Eq. (6) 
is now revised as per Ref. l50 |; th e parameter 9 has been 
determined to be 1039/4620 ^3-) 

The quantity is an undetermined regularization pa- 
rameter that enters the GW flux at 3PN order [8|. An- 
other regularization parameter, uj s , enters the PN expres- 
sions of £ [Eq. i|10|)] and T at 3PN order, and it has been 
determined in the ADM gauge 0, , but not yet in the 
harmonic gauge. However, Eq. © does not depend on 
w s . As in BCV2, we do not include the (partial) spin con- 
tributions to Co at 2.5PN, 3PN, and 3.5PN orders, which 
arise from known 1.5PN and 2PN spin terms of £ and T. 
[To be fully consistent one should know the spin terms of 
£ and T at 2.5PN, 3PN and 3.5PN order.] In Sec- IETD! 
we shall briefly comment on the effect of these terms. We 
ignore also the quadrupole-monopole interaction, which 
we discuss in Sec. IIII El 

The precession equation for the spin is p^HTsf 

Si = o7?( Mw ) 5/3 f 4 + 3 — ) Ljv x Sx , (7) 
2M \ mi ) 

where we have replaced r = r and |Ljv| by their leading- 



order Newtonian expressions in lu, 

r=(j^V 3 , \L n \=^uo = t 1 M^oj- 1 I\ (8) 

The precession of the orbital plane (defined by the normal 
vector Ljv) can be computed following Eqs. (5)-(8) of 
BCV2, and it reads 

i N = ^tt ( 4 + 3— ) Si x L N = Sl L x t N . (9) 

2M \ mi J 

Equations |JS}, J7J), and © describe the adiabatic evo- 
lution of the three variables lu, Si and Ljy- From those 
equations it can be easily deduced that the magnitude 
of the spin, Si — |Si|, and the angle between the spin 
and the orbital angular momentum, ki = Ljv • Si, are 
conserved during the evolution. 

The integration of Eqs. ©, and © should be 
stopped at the point where the adiabatic approximation 
breaks down. This point is usually reached (e. g., for 2PN 
and 3PN orders) when the orbital energy £„pn reaches 
a minimum d£ n ^/dw = (exceptions occur at Newto- 
nian, 1 PN and 2.5 PN orders, as we shall explain in more 
detail in Sec. IIII D|) . We shall call the corresponding orbit 
the Minimum Energy Circular Orbit, or MECO. Up to 
3PN order, and including spin-orbit effects up to 1.5PN 
order, the orbital energy £(u>) reads |6l36j 2 



5 



2 Equation (11) of BCV2 suffers from two misprints: the spin-orbit 
and spin-spin terms should both be divided by M 2 . 



£(o,) = -|(Mu,) 2 /3 



1 - £ + 2) (M.) 2 /3 
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Ljy • Si 
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1 

24 v 
155 



(M^)-^(81-57 ?/ + 77 2 )(Afa;) 4 / 3 
35 



96 



V 



5184 



(Jlfu) 



(10) 



Henceforth, we assume uj s = 0, as computed in Ref. 

B. The precessing convention 

BCV2 introduced a new convention to express the 
gravitational waveform generated by binaries of spin- 
ning compact objects, as computed in the quadrupolar 
approximation; here we review it briefly. At the mass- 
quadrupole leading order, the radiative gravitational field 
emitted by the quasicircular binary motion reads 

*" -§(")<*■ 

where D is the distance between the source and the 
Earth, and Q l J is proportional to the second time deriva- 
tive of the mass-quadrupole moment of the binary, 

Q*i = 2 [A* X j - n l n j ] , (12) 

with n l and A* the unit vectors along the separation vec- 
tor of the binary r, and along the corresponding relative 
velocity v. In general, these vectors can be written as 

n(t)=ei(t) cos$(t) +e 2 (t) sin$(i), (13) 
A(i) = -ei(t) sin$(t) + e 2 (t) cos$(i), (14) 

where &i(t), e2(t), and e^it) = L/y(i) are orthonormal 
vectors, and ei^f) forms a basis for the instantaneous 
orbital plane. 



The adiabatic condition for a sequence of quasi- 
spherical orbits states that n = ui\, but in general $ ^ to. 
The precessing convention introduced by BCV2 is defined 
by imposing that this condition is satisfied (i. e., that 
$ = uj), and it requires that ei j2 (i) precess alongside La? 
as 

e<(i) =n e (t) x e,-(t), i = l,2, (15) 

with 

O e (rj) = L -(0 L -L w )LAr (16) 
[see Eq. © for the definition of fi^] . In this convention, 
the tensor Q Z J can be written as 

Ql j = -2([e+P cos2($ + $ ) + [e x ] 13 sin2($ + $„)) , 

(17) 

with $o an arbitrary initial phase (see below), and 
e + = ei ® ei - e 2 ® e 2 , e x = ei ® e 2 + e 2 ® ei . (18) 



C. The detector response 

The response of a ground-based interferometric detec- 
tor to the GW signal of Eq. (|llf) is given by 



2/iM 

"d7 



[e+P cos2($ + $ ) + [e x ] y sin2($ + *„) [T+] y - F+ + [T x ] y - F x ; (19) 



factor Q: wave generation factor P: detector projection 

the tensors [T +)X ]ij are defined by 

T+=e*®e«-e*®e*, T x = e« ® e« + e« ® e« , (20) 
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after we introduce the radiation frame 

e f = ~ef am<p + e$cosip, (21) 

= — ef cos cos ip — cos sin y> + ef sin , (22) 

ef = +ef sin cos ^ + e£ sin 8 sin <p + ef cos© = N, (23) 

where the detector lies in the direction N with respect to the binary [for the definitions of the angles and tp see 
Fig. 1 of BCV2]. For the antenna patterns F + and F x we have 

F+, x = -[e x <g> e x - e y <g) e y ] ij [T+, x ]y , (24) 
where e x -2/ are the unit vectors along the orthogonal interferometer arms. More explicitly |25|. 

F + = -(1 + cos 2 0) cos20 cos 2-0 - cosQ sin 20 sin 2-0 , (25) 

F x = -(1 + cos 2 0) cos 20 sin 2?/; + cos6> sin 20 cos 2^ (26) 



[for the definitions of the angles 0, and ip, please see 
Fig. 2 of BCV2]. 

Mathematically, we see that the factor P of Eq. (|19fl . 
which is independent of time, collects only terms that de- 
pend on the position and orientation of the detector, and 
that describe the reception of GWs; while factor Q col- 
lects only terms that depend on the dynamical evolution 
of the binary, and that describe the generation of GWs 
(at least if the vectors e 12 .3 arc defined without reference 
to the detector, as we will do soon). Using the language 
of BCV2, in the precessing convention the directional pa- 
rameters 0, (p, 0, 9, and -0 are isolated in factor P, while 
the basic and local parameters of the binary are isolated 
in factor Q. 

Physically, we see that factor Q evolves along three 
different timescales: (i) the orbital period, which sets 
the GW carrier frequency 2$ = 2lu; (ii) the precession 
timescale at which the e +jX change their orientation in 
space, which modulates the GWs; (iii) the radiation- 
reaction timescale, characterized by uj/lu, which drives 
the evolution of frequency. In the adiabatic regime, the 
orbital period is the shortest of the three: so for conve- 
nience we shall define the (leading-order) instantaneous 
G W frequency /gw directly from the instantaneous or- 
bital frequency u>: 

f GW = (2u)/(2n)=u/ir. (27) 

Thus, what parameters are needed to specify Q com- 
pletely? Equation (JSJ for uj(t) can be integrated numer- 
ically, starting from an arbitrary w(0), 3 after we specify 



3 When templates are used in actual GW searches, the initial or- 
bital frequency £00) must be chosen so that most of the sig- 
nal power (i. e., the square of the optimal signal to noise) is 
accumulated at GW frequencies higher than the corresponding 



the basic parameters M, 77, and \ii an d the local param- 
eter K\ = Iin • Si, which is conserved through evolution. 
With the resulting u>(t) we can integrate Eqs. (J7J and 
and then Eq. (|15f) . For this we need initial conditions for 
Si, Ln, and for the e^: without loss of generality, we can 
introduce a (fixed) source frame attached to the binary, 

ef ocSi(0)-[Si(0)-L N (0)]L N (0), 
< = L N (0) xef , ef = L N (0), 
and then take 

ei(0)=ef, e 2 (0) = e^, e 3 (0)=ef. (29) 

[If Si(0) and Ln(0) are parallel, ef can be chosen to lie 
in any direction orthogonal to Ln(0).] The initial orbital 
phase $0 that enters the expression of Q is defined by 

n(0) =ei(0)cos$ + e 2 (0)sin$o, (30) 

while the initial conditions for Si and Ln 7 as expressed 
by their components with respect to the source frame, 
are 

L N (0) = (0,0,1), (31) 
Si(0)= (v/i-k?,0,ki). (32) 



fcw(0) = w(0)/tt, 

y+« h*(f)h(f) [+°° h*(f)h(f) 

For the range of binary masses considered in this paper, and for 
the LIGO-I noise curve, such a /gw(0) should be about 40 Hz. 
Most of the calculations performed in this paper (for instance, 
the convergence tests and the calculation of the mismatch met- 
ric) set instead /gw(0) = 60 Hz to save on computational time; 
experience has proved that the results are quite stable with re- 
spect to this change. 
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FIG. 1: Ending frequency (instantaneous GW frequency at 
the MECO) as a function of n, evaluated from Eq. (1101 at 
2PN order for M = 15Mq, xi = 1| an< l f° r different values of 
Ki. 



BCV2 proposed to use the family of waveforms (detector 
responses) defined by Eqs. ©, 0, ©, GSJ), an d EH as 
a family of physical templates for compact binaries with 
a single spin. Depending on the maximum PN order N 
up to which the terms of Eq. © are retained, we shall 
denote this class of template families STjv- The ST^v 
templates deserve to be called physical because they are 
derived from a physical model, namely the adiabatic PN 
dynamics plus quadrupole GW emission. Each STV tem- 
plate family is indexed by eleven parameters: M, 77, xi 
(basic), Ki (local), 9, <p, 6, <j), i\> (directional), plus the 
initial frequency w(0) (or equivalently, the time to at an 
arbitrary GW frequency), and the initial phase $o- Of 
these, using the distinction between intrinsic and extrin- 
sic parameters introduced in Ref. |26|| 4 and further dis- 
cussed in BCV2, the first four are intrinsic parameters: 
that is, when we search for GWs using STjy templates, we 
need to lay down a discrete template bank along the rele- 
vant ranges of the intrinsic dimensions. The other seven 
are extrinsic parameters: that is, their optimal values 
can be found semialgebraically without generating mul- 
tiple templates along the extrinsic dimensions (another 
way of saying this is that the maximization of the overlap 
over the extrinsic parameters can be incorporated in the 
detection statistic, which then becomes a function only 
of the intrinsic parameters). In Sec. II VI we shall describe 
how this maximization over the extrinsic parameters can 
be achieved in practice. 



D. Comparison between different Post Newtonian 
orders and the choice of mass range 

In this section we investigate the range of masses mi 
and 777.2 for which the PN-expanded evolution equations 
©,0, and © [and therefore the template family (flfjj) ] 
can be considered reliable. As a rule of thumb, we fix 
the largest acceptable value of the total mass by requir- 
ing that the GW ending frequency (in our case, the in- 
stantaneous GW frequency at the MECO) should not 
lie in the frequency band of good detector sensitivity 
for LIGO-I. Considering the results obtained by com- 
paring various nonspinning PN models IUbcvi], and 
considering the variation of the ending frequency when 
spin effects are taken into account [BCV2], we require 
M < 15M Q . In keeping with the focus of this paper 
on binaries with a single significant spin, we also impose 
mijm\ < 0.5, which constrains the spin of the less mas- 
sive body to be relatively small (of course, this condition 
is always satisfied for NS-BH binaries). As a matter of 
fact, population-synthesis calculations j3^ suggest that 
the more massive of the two compact bodies will have the 
larger spin, since usually it will have been formed first, 
and it will have been spun up through accretion from the 
progenitor of its companion. For definitcncss, we assume 
ni\ = 1-3Mq and 777,2 = 7-12M Q ; the corresponding- 
range of 77 is 0.07-0.16. 

In Fig. ^ we plot the GW ending frequency as a func- 
tion of 77, evaluated from Eq. (|10fl at 2PN order for 
M = 15Mq and xi = !• The various curves refer to 
different values of K\. The minimum of the GW ending 
frequency is ~ 300 Hz, and it corresponds to a (12+1)M Q 
binary with spin antialigncd with the orbital angular mo- 
mentum. In Fig. [21 we plot ui/u> 2 , normalized to its lead- 
ing (Newtonian) term 96/5?7(M(jj) 5 / 3 , as a function of 
the instantaneous GW frequency; lu/uj 2 is evaluated from 
Eq. © at different PN orders, for a (10+1.4) M Q binary 
with xi = 1- We see that the effects of the spin-orbit 
interaction (evident for different K\ within the same PN 
order) are comparable to, or even larger than, the effect 
of increasing the PN order. We see also that the differ- 
ent PN curves spread out more and more as we increase 
M and 77. For comparison, in Fig. [3] we show the same 
plot for a (1.4+1.4)M NS-NS binary; note the different 
scale on the vertical axis. In this case the various curves 
remain rather close over the entire frequency band. 

Another procedure (often used in the literature) to 
characterize the effects of spin and PN order on the evolu- 
tion of the GW frequency is to count the number of GW 
cycles accumulated within a certain frequency band: 



A/"gw = 



. dui . 

OJ 



(33) 



4 Note that the concept of extrinsic and intrinsic parameters had 
bee n p resent in the data- ana lysis literature for a long time (see, e. 
g., Sathyaprakash |22| draws the same distinction between 

kinematical and dynamical parameters. 



Here we take cJ m i n = n x 10 Hz and cj max = c^igco = 
(6 3 / 2 7rM) -1 , corresponding to the orbital frequency 
at the innermost stable circular orbit (ISCO) of a 
Schwarzchild black hole with mass M. In Table [I] we 
show A/gw a t increasing PN orders, for (10 + 1.4)M Q , 
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FIG. 2: Plot of e = (Cu /to 2 )/ (96 /5v(M io) 5/3 ) as a function of / G w = w/7T, evaluated from Eq. @ at different PN orders for a 
(10 + 1.4) Mq binary. We do not show the 3.5PN curves, which are very close to the 3PN curves. 
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FIG. 3: Plot of e = (dj /lu 2 )/ (96/5tj(M <o) 5/3 ) as a function of 
/gw = w/7T, evaluated from Eq. © at different PN orders for 
a (1.4 + 1.4) M Q NS-NS binary. We do not show the 2.5PN, 
3PN (9 = 0), and 3.5PN curves, which are very close to the 
2PN curves. Note the change in scale with respect to Fig. 



(12 + 3)M , and (7 + 3)M Q binaries. The contributions 
in parentheses are partial spin terms present at 2.5PN, 
3PN, and 3.5PN orders, and due to known 1.5PN spin 
terms in the orbital energy and luminosity. These terms 
were neglected in Eq. JSJl to be consistent in including 
PN terms of the same order, and we list them here only 
to give their order of magnitude. Unless there are can- 
cellations, the large number of cycles suggests that it is 
worth to compute spin effects up to the 3.5PN order. 

The number of accumulated GW cycles A/qw can be a 
useful diagnostic, but taken alone it provides incomplete 
and sometimes even misleading information. There are 
three reasons for this. First, A/gw is related only to the 
number of orbital cycles of the binary within the orbital 
plane, but it does not reflect the precession of the plane, 
which modulates the detector response in both ampli- 
tude and phase. These modulations are very important 
effects, as witnessed by the fact that neither the stan- 



dard nonspinning-binary templates (which do not have 
builtin modulations) nor the original Apostolatos tem- 
plates (which add only modulations to the phase) can 
reproduce satisfactorily the detector response to the GWs 
emitted by precessing binaries. Second, even if two sig- 
nals have A/"gw that differ by ~ 1 when w max equals the 
GW ending frequency (which apparently represents a to- 
tal loss of coherence, and hence a significant decrease in 
overlap), one can always shift their arrival times to obtain 
higher overlaps. Third, in the context of GW searches the 
differences in A/"gw should be minimized with respect to 
the search parameters, a la fitting factor. 

The Cauchy criterion [271 implies that the se- 
quence STtv converges if and only if, for every k, 
(STjv+fe, STjv) — ► 1 as N — > oo. One requirement of 
this criterion is that (STjv+o.5, STjy) — > 1 as N — > oo, 
and this is what we test in Table ITTI for maximally spin- 
ning and nonspinning (10 + 1.4)M and (12 + 3)M Q bi- 
naries. The overlaps quoted at the beginning of each 
column are maximized over the extrinsic parameters to 
and <&o, but not over the five extrinsic directional param- 
eters if, 0, 9, 4> an( A ip or the intrinsic parameters mi, 
m2, Xi an( i K i- By contrast, we show in parentheses the 
overlaps maximized over all the parameters of the lower- 
order family (i. e., the fitting factors FF for the target 
family STjy+fc as matched by the search family STjy); we 
show in brackets the parameters at which the maximum 
overlaps are achieved. [The overlaps are especially bad 
when IPN and 2.5PN waveforms are used. These two 
orders are rather particular: the flux function T can be 
a decreasing function of u>, and even assume negative val- 
ues (which is obviously not physical) ; correspondingly, to 
can become negative. Furthermore, the MECO criterion 
used to set the ending frequency can also fail, because for 
some configurations the MECO does not exist, or occurs 
after Co has become negative. To avoid these problems, we 
stop the numerical integration of the equations of motion 
when Co decreases to one tenth of its Newtonian value, or 
at a GW frequency of 1 kHz, whichever comes first. For 
comparison, in Table [ri] we show also the overlaps be- 
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(10 + 1.4)M 


(12 + 3)M 


(7 + 3)M Q 


Newtonian 


3577.0 




1522.3 




2283.8 




1PN 


213.1 




114.3 




139.0 




1.5PN 


-181.3 + 




-99.7 + 


55.7ki X i 


-102.3 + 


48.2k iX i 


2PN 


9.8 




6.3 




6.4 




2.5PN 


-20.4 + ( 


21.1kj.Xi) 


-12.7 + ( 


12.1«iXi) 


-10.9 + ( 


9.0«ixi) 


3PN 


2.2 + ( 


-17.0/SiXl + 2.4*1x1) +0.426 


2.2 + ( 


-9.7/dXi + l-2re?x?) +0.406 


2.3 + ( 


-6.6kixi + 0.7rejx?) +0.43 6 


3.5PN 


-1.9 + ( 


6.4kixi) 


-1.3 + ( 


3.8kixi) 


-1.1 + ( 


2.6kixi) 



TABLE I: PN contributions to the number A/"gw of GW cycles accumulated from ui m i n — n x 10 Hz to w max = ojisco = 
l/(6 3,/2 7r M). The contributions in parentheses, "(...)", are partial spin terms present at 2.5PN, 3PN and 3.5PN orders and 
due to known 1.5PN spin terms in the orbital energy and luminosity. 



tween ST2 and ST3, which are much higher than those 
between ST2 and ST2.5, and than those between ST2.5 
and ST 3 .] 

While the nonmaximized overlaps can be very low, the 
FFs are consistently high (note that this requires extend- 
ing the search into the unphysical template region where 
77 > 0.25 and \i > 1); however, the best-fit search param- 
eters can be rather different from the target parameters. 
This suggests that higher-order PN effects can be reab- 
sorbed by a change of parameters, so the ST^r templates 
can be considered rather reliable for the purpose of de- 
tecting GWs from precessing binaries in the mass range 
examined; however, the estimation of binary parameters 
can suffer from systematic errors. In the rest of this pa- 
per we shall describe and analyze a search scheme that 
uses the ST2 template family. 

A more thorough analysis of the differences between 
the various PN orders would be obtained by compar- 
ing the PN-expanded adiabatic model used in this pa- 
per with PN-resummed adiabatic models (d la Pade |27j1 
and nonadiabatic models (a la effective-one-body [I^). A 
similar comparison was carried out for the nonspinning 
case in Refs. [231.BCV1]. Unfortunately, waveforms that 
include precessional effects are not yet available for the 



PN-resummed adiabatic and nonadiabatic models. 



E. The quadrupole— monopole interaction 

In this section we investigate the effect of the 
quadrupole-monopole interaction, which we have so far 
neglected in describing the dynamics of precessing bina- 
ries. It is well known |4G| that the quadrupole moment 
of a compact body in a binary creates a distortion in 
its gravitational field, which affects orbital motion (both 
in the evolution of ui and in the precession of Ljv)j and 
therefore GW emission; the orbital motion, on the other 
hand, exerts a torque on the compact body, changing 
its angular momentum (i. e., it induces a torqued preces- 
sion). Although the lowest-order quadrupole-monopole 
effect is Newtonian, it is smaller than spin-orbit effects 
and of the same order as spin-spin effects. 

When the the spinning body is a black hole, the equa- 
tions for the orbital evolution and GW templates are 
modified as follows to include quadrupole-monopole ef- 
fects. Eq. ijBJl gets the additional term |4l| 



J 



Quad-Mon 



96 



n (M^) 5 /3 J 5 2 "h r (f ^ _ §i)2 _ i 4/3 



2M v ; 



and 



2 A1 M 2 



while the precession equations 0H© become ^lj 



4 + 3^) -3 X i(A^) 1/3 (LjySi) 
mi 



(Ljv x Si) , 



(34) 



(35) 



->JV 



2M 



4 + 3^) -3*i (Mo;) 1 ' 3 (L N ■ Si) 
mi 



furthermore, the orbital energy Hl()|) gets the additional term 



(Si x Ljy) = fl' L x L N ; 



3(Lat-Si) 2 -1 (Mu) 



\4/3 



(36) 



(37) 



r 



last, f2 e is again obtained from Eq. l|16f) . using the mod- ified fl' L in Eq. (|36|l . 
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(N + k,N) 



(ST, 



ki = 0.9 



fc ,STiv> for (10 
ki = 0.5 



+ 1.4)M© binary, xi = 1 
Ki = -0.5 



«i = -0.9 



(1,0) 

(1.5,1) 

(2,1.5) 

(2.5,2) 

(3,2) 

(3,2.5) 

(3.5,3) 



0.1976 (0.7392) [24.5,0.02,0.00,0.00] 

0.2686 (0.7848) [4.53,0.54,0.00,0.00] 

0.4876 (0.99) [9.56,0.14,0.83,0.93] 

0.1587 (0.9578) [10.5,0.13,1.56,0.95] 

0.4395 (0.9848) [11.5,0.10,0.84,0.81] 

0.1268 (0.9758) [12.8,0.08,0.05,0.98] 

0.9614 (0.99) [11.7,0.10,1.00,0.90] 



0.1976 (0.7392) 
0.2696 (0.7008) 
0.5627 ( 0.99) 
0.2011 (0.9887) 
0.5057 (0.9881) 
0.1539 (0.99) 
0.9738 (0.99) 



0.1976 (0.7392) 
0.2065 (0.6040) 
0.6623 (0.99) 
0.2902 (0.9398) 
0.5575 (0.9712) 
0.2520 (0.9744) 
0.9907 (0.99) 



[24.5,0.02,0.00, 
[6.58,0.36,0.00, 
[11.7,0.10,0.97, 
[10.2,0.13,2.00, 
[12.0,0.10,0.92, 
[25.6,0.03,0.35, 
[11.3,0.11,1.02, 



-0.00] 
-0.00] 
-0.50] 
-0.19] 
-0.48] 
-0.21] 
-0.49] 



0.1976 (0.7392) 
0.1800 (0.6255) 



0.7728 
0.3460 
0.6606 
0.2488 
0.9939 



(0.9760) 
(0.99) 
(0.99) 
(0.99) 
(0.99) 



(N + k,N) 



Ki 



0.9 



{ST N+k ,ST N ) for (12 
«i = 0.5 



+ 3)M binary, %i = 1 
Ki = -0.5 



«i 



-0.9 



(1,0) 

(1.5,1) 

(2,1.5) 

(2.5,2) 

(3,2) 

(3,2.5) 

(3.5,3) 



0.2506 (0.7066) [10.5,0 

0.3002 (0.7788) [8.22,0 

0.6379 (0.99) [16.0,0 

0.2039 (0.9397) [15.4,0. 

0.6679 (0.9851) [11.0,0. 

0.1603 (0.99) [18.5,0. 

0.9517 (0.99) [15.2,0 



22,0.00,0.00] 
50,0.00,0.00] 
14,1.14,0.92] 
17,1.95,0.87] 
25,0.72,0.84] 
10,0.05,0.99] 
15,0.84,0.86] 



0.2506 (0.7066) 
0.2597 (0.7381) 
0.7089 (0.99) 
0.2800 (0.9863) 
0.7267 (0.99) 
0.2272 (0.99) 
0.9694 (0.99) 



0.2506 (0.7066) 
0.2124 (0.6934) 
0.8528 (0.99) 
0.4696 (0.9756) 
0.9052 (0.99) 
0.3804 (0.9759) 
0.9932 (0.99) 



[10.5,0.22,0.00, 
[11.6,0.44,0.00, 
[14.2,0.18,1.14, 
[13.5,0.18,1.22, 
[13.1,0.21,1.50, 
[15.8,0.15,0.94, 
[15.3,0.16,1.00, 



-0.00] 
-0.00] 
-0.59] 
-0.51] 
-0.70] 
-0.49] 
-0.49] 



0.2506 (0.7066) 
(0.5427) 
(0.99) 
(0.99) 
(0.99) 
(0.99) 
(0.99) 



0.2017 
0.8620 
0.4219 
0.8868 
0.3060 
0.9900 



(N + k,N) 



(1,0) 

(1.5,1) 

(2,1.5) 

(2.5,2) 

(3,2) 

(3,2.5) 

(3.5,3) 



(STiv+fc, STjv) 
(10 + 1.4) M© 



0.1976 (0.7392) 
0.1721 (0.6427) 
0.7954 (0.9991) 
0.4872 (0.9961) 
0.7471 (0.9970) 
0.4127 (0.9826) 
0.9931 (0.99) 



[24.5,0.02] 
[5.22,0.51] 
[12.7,0.09] 
[6.94,0.25] 
[15.3,0.06] 
[26.5,0.02] 
[11.6,0.11] 



for xi = 

(12 + 3) 



2506 (0.7066) 

2153 (0.6749) 

8924 (0.9981) 

5921 (0.9977) 

8982 (0.9994) 

5282 (0.9783) 

9924 (0.99) 



[10.5,0.22] 
[9.22,0.51] 
[16.2,0.14] 
[8.05,0.48] 
[19.3,0.10] 
[29.0,0.05] 
[15.4,0.15] 



TABLE II: Test of Cauchy convergence of the adiabatic templates STjv at increasing PN orders, for (10 + 1.4)M© and 
(12 + 3)M© binaries, and for maximally spinning BHs (xi = 1, upper and middle panels) and nonspinning BHs (xi = 0, lower 
panel) . The overlaps quoted at the beginning of each column are maximized only with respect to the extrinsic parameters to 
and <E>o- In parentheses, "(...)", we give the overlaps maximized over all the parameters of the lower-order family (i. e., the 
fitting factors FF for the target family STjv+t; as matched by the search family STat, evaluated at the target masses shown); 
the fitting factors are obtained by extending the search into the unphysical template region where r\ > 0.25 and Xi > 1- In 
brackets, "[...]", we show the parameters M, rj,xi,^i (or M, r] when xi — 0) at which the maximum overlap is achieved. The 
detector is perpendicular to the initial orbital plane, and at 3PN order we set 9 — 0; in all cases the integration of the equations 
is started at a GW frequency of 60 Hz. 



The quadrupolc monopole interaction changes the 
number of GW cycles listed in Table |U at 2PN order. 
The additional contributions are 5.2 Xi — 15.5 k\ Xi f° r a 
(10+1.4)M Q binary, 2.5 xf-7.6 k\ X \ for a (12+3)M bi- 
nary, and 1.8x? — 5.4/tf x? f° r a (7 + 3)M Q binary. To es- 
timate more quantitatively the effect of the quadrupolc— 
monopole terms, we evaluate the unmaximized overlaps 
between 2PN templates, computed with and without the 
new terms. The results for (10 + 1.4)M Q binaries are 
summarized in Table IIIII In parentheses we show the 
fitting factors, which are all very high; in brackets we 
show the intrinsic parameters at which the maximum 
overlaps are obtained. We conclude that for the pur- 
pose of GW searches, we can indeed neglect the effects of 
the quadrupole-monopole interaction on the dynamical 
evolution of the binary. 



IV. A NEW PHYSICAL TEMPLATE FAMILY 
FOR NS-BH AND BH-BH PRECESSING 
BINARIES 

In this section we discuss the detection of GWs from 
single-spin precessing binaries using the template family 
first suggested in BCV2, and further discussed in Sec. IIIII 
The proposed detection scheme involves the deployment 
of a discrete template bank along the relevant range of 
the intrinsic parameters M, 77, xi, and K\, and the use 
of a detection statistic that incorporates the maximiza- 
tion of the overlap over all the extrinsic parameters: the 
directional angles 0, ip, 9, (f>, and tp, the time of arrival 
to, and the initial phase $o- In Scc. lIV Al we describe the 
reparametrization of the templates used for the formula- 
tion of the maximized statistic, which is then discussed 
in Sec. II V 51 where we also present an approximated but 
computationally cheaper version. The exact and approx- 
imated statistics are discussed together in Sec. [V] in the 
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view 


«i = 0.9 


(10 + 1.4)Jl 
ki = 0.5 


In* with Vl = 1 

ki = -0.5 


«i = -0.9 


top 
side 

diagonal 


0.4796 (0.99) [10.3,0.13,1-21,0.89] 
0.3503 (0.99) [10.0,0.13,0.77,0.94] 
0.3292 (0.99) [11.2,0.11,0.80,0.94] 


0.9890 (0.99) 
0.8033 (0.99) 
0.6669 (0.99) 


0.1873 (0.99) [11.3,0.11,1.08,-0.48] 
0.8754 (0.99) [11.4,0.11,1.03,-0.39] 
0.4546 (0.99) [11.3,0.11,1.08,-0.49] 


0.7245 (0.9877) 
0.7598 (0.99) 
0.8437 (0.9887) 



TABLE III: Effects of quadrupole-monopole terms, for (10 + 1.4)Af© binaries with maximally spinning BH. At the beginning 
of each column we quote the overlaps between ST2 templates and ST^ M templates that include quadrupole-monopole effects. 
Just as in Tabled these overlaps are maximizing only over the extrinsic parameters to and $o- In parentheses, "(•••)", we show 
the fitting factors for the ST^ M family as matched by the ST2 family; in brackets, "[...]", we show the intrinsic parameters 
at which the fitting factors are achieved. The "view" column describes the position of the detector with respect to the initial 
orbital plane. In all cases the integration of the equations is started at a GW frequency of 60 Hz. 



context of an optimized two-stage detection scheme. 



A. Reparametrization of the waveforms 

We recall from Eqs. |T3 |l -Pfi | that the generic func- 
tional form of our precessing templates is 



h[X 



r[M, 77, xn «i; * , *o5 t]Pij[®, V\ 0, 0, # (38) 



[Please note that for the rest of this paper we shall use 
coupled raised and lowered indices to denote contraction; 
however, the implicit metric is always Euclidian, so co- 
variant and contravariant components are equal. This 
will be true also for the STF components introduced 
later, which are denoted by uppercase roman indices.] 

The factor Q 1 ^ it) (which describes the time-evolving 
dynamics of the precessing binary) is given by 



Q l] = - 



2fi M 



[e+]^cos2($+$ ) + [e x ] 1 



' sin2($+3>o) 
(39) 

where the GW phase $(t) and the GW polarization ten- 
sors e+ iX (t) evolve according to the equations lfl"5|) . l(TB|) 
and (JTSJ. This factor depends on the intrinsic parame- 
ters M, rj, xi, and K\, and on two extrinsic parameters: 
the initial phase <!>0j an d the time of arrival to of the 
waveform, referred to a fiducial GW frequency. We can 
factor out the initial phase $0 by defining 



«tt/2 



Q y (*o = 0), 
Q«($o = tt/4); 



we then have 



QH = Q« cos(2$ ) + Q^ /2 sin(2$ )- 



(40) 
(41) 



(42) 



The factor P^ (which describes the static relative posi- 
tion and orientation of the detector with respect to the 
axes initially defined by the binary) is given by 



the orientation angles 6, <j), and ip, and on the position 
angles and </?, all of them extrinsic parameters. The 
antenna patterns can be rewritten as 



Fx 



= \/F_ 



cos a 
sin a 



(44) 



the factor F = y Ff + F 2 then enters h as an over- 
all multiplicative constant. 5 In what follows we shall be 
considering normalized signals and templates, where F 
drops out, so we set F = 1. We then have 



Pi, 



sin a. 



(45) 



Both Q 1 ^ it) and Py are three-dimensional symmet- 
ric, trace-free (STF) tensors, with five independent com- 
ponents each. Using an orthonormal STF basis M^, 



,5, with iM I ) i jiM J ) ij = 5 



we can convc- 



1 = 1, 

niently express Py and Q* J in terms of their components 
on this basis, 



p 



P J (M J ) 



where 



Q^Q^iM 1 )^, P I = P ij {M I ) 



(46) 



(47) 



In this paper, we shall adopt a particular orthonormal 
basis, 



(MX = 



'-iy 

15 K - y '- 



(M 2 ) 



(M 3 h 



(M 4 



iM% 



22 

ij 



21 

ij 



■ / 47r ^' 

- l VT5 (3: 

V 15 V 1 



15 ^ 
15^' 



(48) 



+yfr 1 ) 



Pij - [F+]ijF+ 



[T x ] ij F x , 



(43) 



where the detector antenna patterns F+. x (0, 4>i "0) an d 
the detector polarization tensors T_|_ lX (©,y>) depend on 



5 In fact, multiple GW detectors are needed to disentangle this 
factor from the distance D to the source. 
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with yf™ defined by 

Xf </V = Y 2m (q), (49) 

where Y 2m (q), m = — 2, . . . ,2 are the usual / = 2 spher- 
ical harmonics, and q is any unit vector. We bring to- 
gether this result with Eqs. I|42l) and (|45|l to write the 
final expression 

h = Pi (Ql cos(2$ ) + Ql /2 sin(2$ )) , (50) 

where 

Pi(G,<p,a) = ([T+(6, </?)]/ cos a + [T x (6, </?)]/ sin a) . 

Henceforth, we shall denote the surviving extrinsic pa- 
rameters collectively as 5 Q = (to, $o> ©> ^P)i an d the 
intrinsic parameters as X 1 = (M, rj, xii K i)' 



B. Maximization of the overlap over the extrinsic 
parameters 



As we have anticipated, it is possible to maximize the 
overlap p = (s, h) semialgebraically over the extrinsic di- 
rectional parameters O, tp, 9, 4>, and ip, without comput- 
ing the full representation of h for each of their configu- 
rations. In addition, it is possible to maximize efficiently 
also over to and $o 5 which are routinely treated as ex- 
trinsic parameters in nonspinning- binary GW searches. 



For a given stretch of detector output s, and for a 
particular set of template intrinsic parameters X 1 = 
(M, i], xi, ki), we denote the fully maximized overlap as 



Ph= = max (s, h(X l , S a )) = max 

H Q to, <&o,@,^, c* 



Pi 



[s, QDto cos 2$ + (s, Qi /2 )t sin 2$ 



PiPj {Qi cos 2$ + sin 2$ , Qi cos 2$ + Qi /2 sin 2$ ) 



(52) 



where the subscript to denotes the dependence of the 
signal-template inner products on the time-of-arrival pa- 
rameter of the templates. In fact, each of these inner 
products can be computed simultaneously for all to with 
a single FFT; in this sense, t is an extrinsic parameter 
0. 

Let us now see how to deal with $ - We start by mak- 
ing an approximation that will be used throughout this 
paper. We notice that the template components PiQ\ 
and PiQ^/2 [Eqs. I|4l)|) and (|41|l ] are nearly orthogonal, 
and have approximately the same signal power, 



over <!>o is found easily: 



Pr 



max 



(s, Qi) cos 2$ + (s, Qi /2 ) sin 2$ 



max 



'PiPjA" 



t ,e,<p, a V PiPjB ij 



^PiP,(Qi,Q{) 



max p$ 



where we have defined the two matrices 



A IJ = 



D 



u _ 



(55) 



(s,Qi)t (s,Qi)t + (s,Qi /2 )t {s,Qi /2 )t , 



(56) 



(PjQlPjQt/z) 
(PiQlPjQi) 



PiQU,PjQl 



/2) 



(53) 
(54) 



this is accurate as long as the timescales for the radiation- 
reaction-induced evolution of frequency and for the 
precession-induced evolution of phase and amplitude 
modulations are both much longer than the orbital pe- 
riod. More precisely, Eqs. i|53|) and (|54|l are valid up to 
the leading-order stationary-phase approximation. Un- 
der this hypothesis Eq. (|52|l simplifies, and its maximum 



which arc functions only of the intrinsic parameters (and, 
for A IJ , of to). We have tested the approximations i|53|) 
and (|54|l by comparing the maximized overlaps obtained 
from Eq. I|55|) with the results of full numerical maximiza- 
tion without approximations; both the values and the lo- 
cations of the maxima agree to one part in a thousand, 
even for systems with substantial amplitude and phase 
modulations, where the approximations are expected to 
be least accurate. 

Although Eq. I|55|) looks innocent enough, the maxi- 
mization of p$ is not a trivial operation. The five com- 
ponents of Pi in Eq. (|55|l are not all independent, but 
they are specific functions of only three parameters, 0, 
(p., and a [see the discussion leading to Eqs. I|45|) and 
Q51JI-] We can therefore think of ps a as the result of 
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maximizing p$ with respect to the five-dimensional vec- 
tor Pj, constrained to the three-dimensional physical sub- 
manifold P/ (0, <p, a). We shall then refer to p^<* as the 
constrained maximized overlap. 

What is the nature of the constraint surface? We 
can easily find the two constraint equations that define 
it. First, we notice from Eqs. 1521) and (|55|l that the 
magnitude of the vector Pi does not affect the overlap: 
so we may rescale Pi and set one of the constraints as 
PjP 1 = 1; even better, we may require that the denom- 
inator of Eq. (J5SJ be unity, PjPjB ij = 1. Second, we 
remember that pj [Eq. I|43|)] is the polarization tensor 
for a plane GW propagating along the direction vector 



N l = (sinG cos <p, sinG simp, cos0). 



(57) 



Because GWs are transverse, Pjj must admit N l as an 
eigenvector with null eigenvalue; it follows that 



detP,, 



0. 



(58) 



This equation can be turned into the second constraint 
for the Pj [see Eq. QKty of App-EJ. 

Armed with the two constraint equations, we can re- 
formulate our maximization problem using the method 
of Lagrangian multipliers [Eq. (|A7|) in App. El- How- 
ever, the resulting system of cubic algebraic equations 
does not appear to have closed- form analytic solutions. 
In App. ^ we develop an iterative algebraic procedure 
to solve the system, obtaining the constrained maximum 
and the corresponding Pj. In practice, we have found it 
operationally more robust to use a closed-form expression 
for the partial maximum over $o and a (which seems to 
be the farthest we can go analytically), and then feed it 
into a numerical maximum-finding routine (such as the 
well-known amoeba 43]) to explore the (0, tp) sphere, re- 
peating this procedure for all to to obtain the full maxi- 
mum. 

To maximize over a, we use Eq. (|51|l to factor out 
the dependence of the Pj on a, and write 



' PiPjA ij 
PiPjB" 



I uA a u T 
uB„u T ' 



(59) 



where u is the two-dimensional row vector (cos a, sin a), 
and where A Q and B Q are the two-by-two matrices 



A 

B,, 



= A 



u 



B 



u 



[T+]i[T+]j 
[T+]/[T x ]. 7 

[T+]/[T + ]. 7 
[T + ]/[T x ], 7 



[T> 



]/[Tx]j 
]/[Tx]./ 

[T+MTx]j 
[Tx]/[T x ]j 



(60) 
(61) 



in these definitions we sum over the indices I and J. The 
matrices A Q and B Q are implicitly functions of the angles 
and ip through the polarization tensors T + and T x . It 
is straightforward to maximize Eq. (|59|) over a, yielding 6 



max 

toi©,V 



maxeigv [A a B Q 1 ] 



max p<s> 0t u 



(62) 



The overlap p$ , Q is essentially equivalent to the T statis- 
tic used in the search of GWs from pulsars |44j . 

The last step in obtaining p-c, is to maximize p$ 0jCt nu- 
merically over the (0, (p) sphere, repeating this procedure 
for all to to obtain the full maximum. Now, to enters Eq. 
l|62Jl only through the ten signal-template inner prod- 
ucts (siQott/2) con tained in A Q , and each such prod- 
uct can be computed for all to with a single FFT. Even 
then, the semialgebraic maximization procedure outlined 
above can still be very computationally expensive if the 
search over and ip has to be performed for each individ- 
ual to. We have been able to reduce computational costs 
further by identifying a rapidly computed, fully algebraic 
statistic p~ a that approximates /?~a from above. We then 
economize by performing the semialgebraic maximiza- 
tion procedure only for the values of to for which pL a 
rises above a certain threshold; furthermore, the loca- 
tion of the approximated maximum provides good initial 
guesses for and ip, needed to kickstart their numerical 
maximization. 

Quite simply, our fast approximation consists in ne- 
glecting the functional dependence of the Pj on the di- 
rectional parameters, computing instead the maximum 
of p$ [Eq. H55fl ] as if the five Pi were free parameters. 
Using the method of Lagrangian multipliers outlined in 
the beginning of App. El [Eqs. (JEMS^], we § et 



with 



(A IJ - \B 1J )Pj = 



max - 

Pi 



IJ\ 



I PiPjA" 
PiPjB ij 



= ^/maxeigv [AB- 1 ] , (63) 



A = maxeigv[AB _1 ] . (64) 



Here the prime stands for unconstrained maximization 
over Pf. We shall henceforth refer to p' Sa as the uncon- 
strained maximum. 

Note that the value of the Pj at the unconstrained 
maximum will not in general correspond to a physical 
set of directional parameters, so P,j will not admit any 
direction vector N l [Eq. I|57|) ] as a null eigenvector. How- 
ever, we can still get approximate values of and ip by 
using instead the eigenvector of Py with the smallest 
eigenvalue (in absolute value). 



6 Just as it happens for the Pj, the magnitude of u does not affect 



the value of Eq. 1591 . so the maximization can be carried out 
cquivalcntly over all the vectors u that satisfy uB tt u T = 1. We 
can then use a Lagrangian-multiplier method to find the maxi- 
mum, Eq. 1621 . and the corresponding u, in a manner similar to 
the procedure used in App. [X] 
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V. DESCRIPTION AND TEST OF A 
TWO-STAGE SEARCH SCHEME 



Numerical comparison of constrained and 
unconstrained maximized overlaps 



In Sec. IIVI we have described a robust computational 
procedure to find the maximum overlap p-a (which is 
maximized over the extrinsic parameters to, and Pj, 
where the allowed values of the Pj are constrained by 
their functional dependence on the directional angles). 
We have also established a convenient analytic approx- 
imation for , the unconstrained maximized overlap 
p S a (which is maximized over the extrinsic parameters 
<f>0i ^Oj and Pj, but where the Pj are treated as five in- 
dependent and unconstrained coefficients). Because the 
unconstrained maximization has access to a larger set of 
P/, it is clear that p' Sa > p-c Still, at least when the 
target signal s is very close to the template h(Xi), we 
expect p'-a to be a very good approximation for p-c, . 

A quick look at the procedures outlined in Sec. IIVI 
shows that, for the filtering of experimental data against 
a discrete bank of templates {h(XJ k ^)}, the computation 
of p' 3a is going to be much faster than the computation 
of ps." ■ Under these conditions, it makes sense to im- 
plement a two-stage search scheme where the discrete 
bank {h(X^)} is first reduced by selecting the templates 
that have high p' Sa against the experimental data; at this 
stage we identify also the promising times of arrival to- 
The exact ps° is computed only for these first-stage trig- 
gers, and compared with the detection threshold p* to 
identify detection candidates (one would use the same 
threshold p* in the first stage to guarantee that all the 
detection candidates will make it into the second stage). 7 

To prove the viability of such a search scheme, we shall 
first establish that pta is a good approximation for pz<* 
for target signals and templates computed using the adi- 
abatic model of Sec. IIIII We will take slightly displaced 
intrinsic parameters for target signals and templates, to 
reproduce the experimental situation where we are trying 
to detect a signal of arbitrary physical parameters with 
the closest template belonging to a discrete bank. This 
first test is described in Sec. IV Al We shall then study the 
false-alarm statistics of p-a and p' Sa , and we shall show 
that, for a given detection threshold, the number of first- 
stage triggers caused by pure noise is only a few times 
larger than the number of bona fide second-stage false 
alarms. Such a condition is necessary because the two- 
stage detection scheme is computationally efficient only 
if few templates need ever be examined in the expensive 
second stage. The false-alarm statistics (in Gaussian sta- 
tionary noise) are obtained in Sec. IV Bl and the second 
test is described in Sec. IV Cl 



In this section we describe a set of Monte Carlo runs 
designed to test how well pt Q can approximate pz<* , for 
the target signals and templates computed using the adi- 
abatic model of Sec. IIIII for typical signal parameters, 
and for signal-template parameter displacements char- 
acteristic of an actual search. 

We choose target signals with 20 different sets of in- 
trinsic parameters given by 



(mi,m 2 ,xi,Ki) 



(1O,1.4)M 
(7,3)M 



0.5 
1 



7 This is not a conventional hierarchical scheme, at least not in the 
sense that there is a tradeoff between performance and accuracy. 



-0.9 
-0.5 

\ < o.o 

0.5 
0.9 

(65) 

For each set of target-signal intrinsic parameters, 
we choose 100 random sets of extrinsic parameters 
(0, <p, a, $o)) where the combination (0, ip) is distributed 
uniformly on the solid angle, and where a and <E>o are 
distributed uniformly in the [0, 2n] interval. The target 
signals are normalized, so the allowed range for p~a and 
pUis[0,l]. 

For each target signal, we test 50 (normalized) 
templates displaced in the intrinsic-parameter space 
{M, T}, xi, Ki) [the optimal extrinsic parameters will be 
determined by the optimization of p~ Q and pg Q , so we 
do not need to set them]. The direction of the displace- 
ment is chosen randomly in the [M, 77, X\i K i) space. For 
simplicity, the magnitude of the displacement is chosen 
so that, for each set of target-signal intrinsic parameters 
and for the first set of target-signal extrinsic parameters, 
the overlap p'~ a is about 0.95; the magnitude is then kept 
fixed for the other 99 extrinsic-parameter sets, so p^ Q can 
be very different in those cases. 

Figure 0] shows the ratio p'~ a / ps° , for each pair [20 x 
50 in total] of target and template intrinsic-parameter 
points, averaged over the 100 target extrinsic-parameter 
points, as a function of the averaged ps." ■ The p' Sa get 
closer to the ps a as the latter get higher; most important, 
the difference is within ~ 2% when p^a > 0.95, which 
one would almost certainly want to achieve in an actual 
search for signals. We conclude that p' Sa can indeed be 
used as an approximation for ps<* in the first stage of 
a two-stage search. The second stage is still necessary, 
because the false-alarm statistics are worse for the un- 
constrained maximized overlap (where more degrees of 
freedom are available) than for its constrained version. 
We will come back to this in the next two sections. 

It is also interesting to compare the set of extrinsic pa- 
rameters of the target signal with the set of extrinsic pa- 
rameters that maximize p~a , as characterized by the cor- 
responding source direction vectors, N truc and N max re- 
spectively. Figure[S]shows the inner product N tr ue'N max , 
averaged over the 100 target extrinsic-parameter points, 
as a function of the averaged p~ Q . The difference between 
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FIG. 4: Ratio between the unconstrained (p=a) and constrained (p~°) maximized overlaps, as a function of p~a. Each point 
corresponds to one out of 20 x 50 sets of intrinsic parameters for target signal and template, and is averaged over 100 sets of 
extrinsic parameters for the target signal. The error bars show the standard deviations of the sample means (the standard 
deviations of the samples themselves will be 10 times larger, since we sample 100 sets of extrinsic parameters). The two 
panels show results separately for (10 + 1.4) Mq (left) and (7 + 3)Mq target systems (right). The few points scattered toward 
higher ratios and lower ps° are obtained when the first set of extrinsic parameters happens to yield a high p' sa that is not 
representative of most other values of the extrinsic parameters; then the magnitude of the intrinsic-parameter deviation is set 
too high, and the comparison between p' sa and p=a is done at low p=» , where the unconstrained maximized overlap is a poor 
approximation for its constrained version. 




FIG. 5: Inner product between target-signal source direction Ntrue and p H a -maximizing source direction N max , as a function 
of ph° ■ Each point corresponds to one out of 20 x 50 sets of intrinsic parameters for target signal and template, and is averaged 
over 100 sets of extrinsic parameters for the target signal. Standard deviations of the sample means are shown as error bars, 
as in Figure[I] The two panels show separately (10 + 1.4) Mo target systems (left) and (7 + 3)Mq target systems (right). 



the vectors can be very large, even when p=:a > 0.95: 
this happens because the intrinsic-parameter displace- 
ment between target signal and template can be compen- 
sated by a change in the extrinsic parameters of template 
(in other words, the effects of the intrinsic and extrinsic 
parameters on the waveforms are highly correlated). 



B. False-alarm statistics for the constrained and 
unconstrained maximized overlaps 

In this section we derive and compare the false-alarm 
statistics of ps" and p' Sa ■ Our purpose is to estimate the 
number of additional triggers that are caused by replac- 



ing the detection statistic by the first-stage statistic 
P^q . Our two-stage detection scheme, which employs the 
rapidly computed p' sa to choose candidates for the more 
computationally expensive ps° , will be viable only if the 
number of those candidates is small enough. 

By definition, a false alarm happens when, with in- 
terferometer output consisting of pure noise, the detec- 
tion statistic computed for a given template happens to 
rise above the detection threshold. Although the detec- 
tion statistics ps." and p'~ a include maximization over the 
time of arrival to , we find it convenient to exclude to from 
this computation, and to include it later when we evalu- 
ate the total false-alarm probability for all the templates 
in the bank. 
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Recall that p^ [Eq. (JSSJ] and p'^ a [Eq. (jHSJ] are func- 
tions of the matrices A and B, which contain the inner 
products (s, Q0.Tr/2) and ( ( 3o,7r/2' respectively. In 

this case the signal s is a realization of the noise, n. Wc 
combine the vectors Qq and <2^y 2 together as Q x with 
X = 1, . . . , 10; under the assumption of Gaussian station- 
ary noise, Y x = (n, Q x ) is a ten-dimensional Gaussian 
random vector with zero mean and covariance matrix 1.351 



the marginal probability density 



G IJ = (n,Q x )(n,QJ) = (Q X ,Q J ) 



(66) 



The covariance matrix C XJ specifies completely the sta- 
tistical properties of the random vector Y , and it is 
a function only of B, and therefore only of the intrin- 
sic parameters of the template. We can also combine 
P/cos2$o and Pi sin2<E>o together as P x , and then write 
the maximized overlaps p~. a and p' sa as 



Px(n,Q J ) 



P1Y 1 



(67) 



where maximization if performed over the appropriate 
range of the P x . In the rest of this section we shall use 
the shorthand p to denote both p^a and p' Sa . 

Equation l|67|l is very general: it describes p= a and p'~ a , 
but it can also incorporate other maximization ranges 
over the P x , and it can even treat different template 
families. In fact, the maximized detection statistic for 
the (^/'oV'3/2'B)6 DTF of BCV2 can be put into the same 
form, with Pj = aj, for X — 1, . . . , 6, and with com- 
pletely unconstrained maximization. 

We can now generate a distribution of the detection 
statistic p for a given set of intrinsic parameters by gen- 
erating a distribution of the Gaussian random vector Y x , 
and then computing p from Eq. (|67|l . The first step is 
performed easily by starting from ten independent Gaus- 
sian random variables Z x of zero mean and unit variance, 
and then setting Y x = [y/C] XJ Zj. 8 Thus, there is no 
need to generate actual realizations of the noise as time 
series, and no need to compute the inner products (n, Q x ) 
explicitly. 

The statistics p [Eq. 1|67|) ] are homogeneous with re- 
spect to the vector Z x : that is, if we define Z 1 = rZ x 
( where r = y Z X Z X and Z X Z X = 1) we have 



p[Y J -(Z x )}=rp[Y x (Z x )}=rp 1 (n) 



(68) 



here f2 represents the direction of Z x in its ten- 
dimensional Euclidian space. The random variable r has 



The square root of the matrix \\fG\ L ^ can be defined, for in- 
stance, by \/C\/C T = C, and it can always be found because 
the covariance matrix C 1 ^ is positive definite. It follows that 
[v / C] I£ [v / Cp- M MX7 



yiyj 

C x 3 , as required 



Pr(r) 



r"- 1 exp(-r 2 /2) 
2"/ 2 - 1 r(i//2) 



(69) 



where the direction f2 is distributed uniformly over a ten- 
sphere. [For the rest of this section we shall write equa- 
tions in the general ^-dimensional case; the special case 
of our template family is recovered by setting v = 10.] 
The random variables r and f2 [and therefore pi(Q)] are 
statistically independent, so the cumulative distribution 
function for the statistic p is given by the integral 



P(p < p*) = dn 



p'/piip.) 



„ (pVpi(Q)) 

2 ' 2 



p r (r) dr J J d£l 

dn/ldn, (70) 



where T[fc,z] = J z + °° t k ~ 1 e~ t dt is the incomplete gamma 
function. 

The false-alarm probability for a single set of intrin- 
sic parameters and for a single time of arrival is then 
1 — P(p < p*). The final integral over the i^-dimensional 
solid angle can be performed by Monte Carlo integration, 
averaging the integrand over randomly chosen directions 
f). Each sample of the integrand is obtained by generat- 
ing a normalized Z x (that is, a direction f2), obtaining 
the corresponding Y x , computing pi(fl) from Eq. 1)67(1 . 
and finally plugging pi(£l) into the T function. 

Equation l|70|l shows that if we set pi(£l) = 1, the ran- 
dom variable p follows the \{v) distribution; this is obvi- 
ous because in that case p — r = ^JZ x Zx [see Eq. £3)], 
where the Z x are v independent Gaussian random vari- 
ables. In fact, jOi(fi) can be written as 



R X Z 2 



where R x = [VC] XJ P J ; 

. 

which shows that pi (11) = 1 uniformly for every Q if and 
only if the range of maximization for P x is the entire 
i/-dimensional linear space generated by the basis {Q x }; 
however, once we start using the entire linear space, the 
particular basis used to generate it ceases to be impor- 
tant, so the covariance matrix C x ^ drops out of the equa- 
tions for the false- alarm probabilities. That is the case, 
for instance, for the {tpo^^^B)^ DTF [see Sec. V B of 
BCV2], whose false-alarm probability is described by the 
X(i/=6) distribution. For our template family v = 10, 
but both p~,a and p' Sa have very restrictive maximiza- 
tion ranges for P x (because Pz=i,...,5 and Pi=6,...,io are 
strongly connected): so both p-~a and p'~ a will have much 
lower false-alarm probability, for the same threshold p* , 
than suggested by the X(u=io) distribution. In fact, in 
the next section we shall see that the effective v for the 
detection statistic p' 3a is about 6; while the effective v 
for p~a is even lower. 
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C. Numerical investigation of false-alarm statistics 

The total false-alarm probability for the filtering of 
experimental data by a template bank over a time T is 

P tot (p > p* ) = 1 - [P(p < p* )] Ar ~ Ar «- (72) 

(see for instance BCV1), where the exponent 
A/"shapesA/timcs is an estimate of the number of ef- 
fective independent statistical tests. The number of 
independent signal shapes A4hapos is related to (and 
smaller than) the number of templates in the bank; 9 the 
number of independent times of arrival .A/times is roughly 
T/6t , where dto is the mismatch in the time of arrival 
needed for two nearby templates to have, on average, 
very small overlap. In our tests we set A/" s hapcs = 10 6 
and A/times = 3 x 10 10 (or equivalently St ~ 1 ms), as 
suggested by the results of Sec. El for template counts 
and for the full mismatch metric; in fact, both numbers 
represent rather conservative choices. 

We compute single-test false-alarm probabilities from 
Eq. (|70|l . averaging the integrand over 10 5 randomly cho- 
sen values of fl to perform the integration over Q, a la 
Monte Carlo. Our convergence tests indicate that this 
many samples are enough to obtain the required preci- 
sion. 10 In Fig. we show the thresholds p* required 
to achieve a total false-alarm rate of 10 _3 /year; the fig- 
ure suggests that a threshold close to 10 is adequate. 
The thresholds are only marginally higher for the un- 
constrained statistic, so the number of first-stage false 
alarms that are dismissed in the second stage is lim- 
ited. We show also the threshold required to achieve 
the same false-alarm rate with the (ipo ^3/2-8)6 DTF of 
BCV2: this threshold is very close to the values found 
for p' Sa , indicating that p' Ea has roughly six effective de- 
grees of freedom (as it seems reasonable from counting 
the five P 1 plus $ )- The BCV2 threshold is consistently 
higher than the ps^ threshold for the same single-test 
false-alarm rate; this suggests that the detection scheme 
discussed in this paper is less wasteful (with respect to 
the available signal power) than the BCV2 scheme, as- 
suming of course that the number of templates used in 
the two banks is similar. 

In Fig. we show the ratio between the single-test 
false-alarm probabilities for p-~a and p'~ a : for a common 
threshold around 10, we can expect about five times more 
false alarms using p' 3c , than using p~<* , for most values of 
the intrinsic parameters (for some of them, this number 
could be as high as ~ 15). These results corroborate 



It is given, very roughly, by the number of templates when the 
minimum match is set to 0. See for instance BCV1. 
In fact, the average is dominated by the samples that yield 
the larger values of p±(Q), since the V function amplifies small 
changes in its argument. So the number of samples used in the 
Monte Carlo integration needs to be such that enough large pi (£2) 
do come up. 
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FIG. 6: Detection thresholds p* for a false-alarm rate of 
KTVyear, 

using the constrained statistic p~° (circles), the 
approximated (unconstrained) statistic p' B a (diamonds), and 
the detection statistic for the (ipotp 3 / 2 B)6 DTF from BCV2 
(dashed line) . The four windows correspond to the masses and 
Xi shown; the points in each window correspond to m = 0.9, 
0.5, 0, -0.5, -0.9. 



our suggestion of using p' Ea in the first-stage of a two- 
stage detection scheme, to weed out most of the detection 
candidates before computing the more computationally 
expensive p^a . 



VI. TEMPLATE COUNTING AND 
PLACEMENT 

The last aspect to examine before we can recommend 
the template family of Sec. IIVI for actual use with the 
two-stage search scheme of Sec.[V]is the total number of 
templates that are needed in practice. As mentioned in 
Sec.^ the template-bank size and geometry required to 
achieve a certain minimum match can be studied using 
the mismatch metric prl l28l l33| . which describes, to 
quadratic order, the degrading overlap between nearby 
elements in a template bank: 

1- (h(\ A ),h(\ A + A\ A )) 

= S[X A , X A + AX A ] = g BC AX B AX c , (73) 

where 6 denotes the mismatch, and where 

i d 2 (h(x A ),h(x A + Ax A )) 

9BC 2 d(AX B )d(AX c ) ' [ ' 

No zeroth- or first-order terms are needed in the expan- 
sion (|73|) . because the overlap has a maximum of 1 (for 
normalized templates) at AX A = 0. The metric is posi- 
tive definite, because 8 > 0. Note that, according to this 
definition, the mismatch 5 is the square of the metric 
distance between X A and X A + AX A . It is also half the 
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FIG. 7: Ratio (1 — < p*])/(l — -P[p~° < P*\) between single-test false-alarm probabilities for the unconstrained and con- 

strained detection statistics, as a function of threshold p* . The two panels represent systems with masses equal to (10+1. 4)M© 
(left) and to (7+3)Mq (right). The five curves in each plot correspond to different Ki- 



square of the inner-product distance y (Ah, Ah), where 

Ah = h(\ A ) - h(X A + AA A ). n 

Ideally, for a given continuous template family, one 
could find a reparametrization in which the metric is a 
Kronecker delta, and then lay down a template bank as a 
uniform hypercubic lattice in these coordinates, with the 
appropriate density to yield the required MM. For a hy- 
percubic lattice in n dimensions, 12 the (metric) side 51 of 
the lattice cell is given by the relation 1 — MM = n(Sl/2) 2 
[2lll2^ |: we then get the total number of templates in the 
bank by dividing the total (metric) volume of parameter 
space by the volume of each cell: 



templates 



V\detg BC \d n X A 2^/(1 - MM)/n 



(75) 

In practice, this expression will usually underestimate 
the total number of templates, for two reasons: first, for 
more than two dimensions it is usually impossible to find 
coordinates where the metric is diagonalized everywhere 
at once; second, the fact that the actual parameter space 
is bounded will also introduce corrections to Eq. (|75|1 . 
[The presence of null parameter directions, discussed in 
Sec. I VI lH can also be seen as an extreme case of bound- 
ary effects.] 

As we showed in Sees. IIVI and the overlap of the 
detector output with one of the STjv templates can be 
maximized automatically over all the extrinsic parame- 
ters E a ; it follows that a discrete template bank will need 
to extend only along the four intrinsic parameters X 1 . So 
the estimate <f75|) for the number of templates should be 



11 So gBC is truly a metric, except for a factor of 1/2. 

12 For specific dimensionalities, other regular packings might be 
more efficient: for instance, in two dimensions a lattice of equilat- 
eral triangles requires fewer templates than a lattice of squares. 



computed on the projected metric gf, roj that satisfies 



l-ps« h(X\ E a ), h(X i + AX ! ) 
= 1 - mnx(h(X\~ a ),h(X l + AX\E' a )) (76) 



? P roj AI l AP. 



Note that <?f 7 roj is still a function of all the parameters. 



In Sec. IVI Al we compute gfj 0} from the full metric gsc\ 

we then proceed to construct an average metric, gf™ : ', 
which is connected closely to detection rates, and does 
not depend on the extrinsic parameters. 

In fact, it turns out that not all four intrinsic param- 
eters are needed to set up a template bank that achieves 
a reasonable MM: we can do almost as well by replacing 
a 4-D bank with a 3-D bank where (for instance) we set 
K\ =0. As a geometrical counterpart to this fact, the 
projected metric must allow a quasinull direction: that 
is, it must be possible to move along a certain direction 
in parameter space while accumulating almost no mis- 
match. The correct template counting for the 3-D bank 
is then described by a reduced metric, which we discuss 
in Sec. IVIES! Finally, we give our results for the total 
number of templates in Sec. IVI CI 



A. Computation of the full, projected and average 
metric 



According to Eq. Ij74|) . the full metric gsc can be com- 
puted numerically by fitting the quadratic decrease of the 
overlap (h(\ A ),h(\ A + A\ A )} around AX A = 0. It is also 
possible to rewrite gsc i n terms of first-order derivatives 
of the waveforms: since (h(\ A ), h(X A )) — 1 for all \ A , 

9 <M> = 2(/>,|^)=0 (77) 
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[in this equation and in the following, we omit the para- 
metric dependence h(X A ) for ease of notation]; taking one 
more derivative with respect to X c , we get 



dh dh 



dX c ' dX B 



which implies [by Eq. I|74|l ] 



h, 



d 2 h \ 
dx c d\ B J 



o. 



1 / dh dh \ 
9BC ~ 2\9AS' d>P/' 



(78) 



(79) 



The inner product in the right-hand side of Eq. I|79|) ex- 
presses the Fisher information matrix for the normalized 
waveforms h(X A ) (see for instance Ref. pjjf'l: for nonnor- 
malized waveforms h(X A ) we can write 



9BC 



dh dh \ 



2(h,h) \dX B ' dX c / 

1 / dh \ I dh_\ 
2(h,h)i \dX B ' / \ 1 dX c I 



(80) 



It is much easier to compute the mismatch metric 
from Eq. (|79|) rather than from Eq. Ij74(l . for two rea- 
sons. First, we know the analytic dependence of the 
templates on all the extrinsic parameters (except to), 
so we can compute the derivatives dh/dE a analyti- 
cally (the derivative with respect to to can be handled 
by means of the Fourier-transform time-shift property 
F[h(t + t )] = .F[/i(t)]exp[27u/to]). Second, although 
the derivatives dh/dX l have to be computed numeri- 
cally with finite-difference expressions such as [h(X i + 
AX l ,E a ) - h(X l ,E a )}/AX\ this is still easier than fit- 
ting the second-order derivatives of the mismatch numer- 
ically. 13 

To obtain the projected metric gfj°\ we rewrite the 
mismatch S(X A ,X A + AX A ) by separating intrinsic and 
extrinsic parameters, 



5{X\E a ;X i + AX\E' 
= ( AX' AE a ) 



fA» Q ) 

Gij 

G a j "fafi 



AX j 
AS 13 



(81) 



13 We have found that we can obtain a satisfactory precision for 
our metrics by taking several cautions: (i) reducing the pa- 
rameter displacement AX 1 along a sequence ( fc )AP until the 
norm (WA/i, - ^"^Ah, ^Ah - V'-^Ah) of the feth correc- 
tion becomes smaller than a certain tolerance, where ^Ah = 
\h(X* + WAX',H a ) - h(X l ,S a )] /(*)AX ! is the fcth approxi- 
mation to the numerical derivative; (ii) employing higher-order 
finite-difference expressions; (iii) aligning both the starting and 
ending times of the waveforms h(X',S a ) and h{X % + AJ*,S") 
by suitably modifying their lengths [by shifting the two wave- 
forms in time, and by truncating or extending h{X % + AX',2°) 
at its starting point]. 



here we have split the full metric gsc into four sec- 
tions corresponding to intrinsic-intrinsic {Gij), extrinsic- 



extrinsic ("/ a p), and mixed (C Q 



Cj a ) components. 



Maximizing the overlap over the extrinsic parameters is 
then equivalent to minimizing Eq. I|81() over the AS" for 
a given AX 1 , which is achieved when 



-YapASP + CajAX* =0. 

while the resulting mismatch is 

mui6(X i ,E a ;X i + AX\E a + AE a ) -- 



(82) 



A2 



[G^-CUl'^Cp,] AX l AX? 



— r,P ro J 



AX'AXi 



(83) 



Here (7 1 ) a @ is the matrix inverse of 7 a /3- For each point 
(X 1 , S a ) in the full parameter space, the projected metric 
gfj 0i describes a set of concentric ellipsoids of constant 
Phq in the intrinsic-parameter subspace. We emphasize 
that the projected metric has tensor indices correspond- 
ing to the intrinsic parameters, but it is a function of 
both the intrinsic and the extrinsic parameters, and so 
are the constant-psc, ellipsoids. 

Therefore, to build a template bank that covers all the 
signals (for all X 1 and S Q ) with a guaranteed MM, we 
must use the projected metric at each X 1 to construct 
the constant-mismatch ellipsoids for all possible S Q , and 
then take the intersection of these ellipsoids to determine 
the size of the unit template-bank cell. This is a minimax 
prescription |27j | , because we are maximizing the overlap 
over the extrinsic parameters of the templates, and then 
setting the template-bank spacing according to the least 
favorable extrinsic parameters of the signal. In general, 
the intersection of constant-mismatch ellipsoids is not an 
ellipsoid, even in the limit 5 — > 0, so it is impossible 
to find a single intrinsic-parameter metric that can be 
used to enforce the minimax prescription. There is an 
exception: the projected metric is not a function of t Q 
or $0i 14 so it can be used directly to lay down banks of 
nonspinning-binary templates |2a , |28| , for which to and 
$0 are the only extrinsic parameters. 

Returning to the generic case, we can still use the pro- 
jected metric to guide the placement of a template bank 
if we relax the minimax prescription and request that 
the minimum match be guaranteed on the average for 
a distribution of signal extrinsic parameters. It turns 
out that this average-mismatch prescription is closely re- 
lated to the expected detection rates. Let us see how. 
The matched-filtering detection rate for a signal s 
S A x h(X l , E a ), where SA = (s, s) 1 / 2 is the signal ampli- 
tude at a fiducial luminosity distance, is proportional to 



The overlap {h(X i , E a ), h(X'\ S'")) depends only on to - 
and <I>o — <t> : we have h(f) ~ exp{27ri/to + i<&o} and h'(f) 



e x p{2nift' + i<S>' }, so h'(f)h*(f) ■ 
<£o)}- 



exp{27ri/(t -t ) + i($ - 
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SA^ysAicar], where /i ncar = h{X l + AX\E' a ) is the 
closest template in the bank, and where we assume that 
sources are uniformly distributed throughout the volume 
accessible to the detector (see, for instance, BCV1). The 
minimax prescription is given by 

PE 4§, h ncar ] ~ 1 - 9%°\X\ Z a )AX l AXi > MM (84) 

for all S", which ensures that the detection rate is re- 
duced at most by a factor MM 3 for every combination of 
signal extrinsic and intrinsic parameters. 

Averaging over a uniform distribution of signal extrin- 
sic parameters, 15 we get a detection rate proportional to 



/ 



gE q SA 3 ( 1 - ,,!';'• A.Y'A.Y-' 



SA 3 



dE a SA- i g, 



*3 



AX 1 AX 3 



SA 3 ( 1 - qff oi AX l AXn 



(85) 



where SA 3 = J Gfc,"SA 3 , and where the average metric 
gfj° 3 , now a function only of X 1 , is defined as 



gfP= I dS Q SA 3 5 r J /SA 3 



(86) 



[To derive Eq. (|S3|) we assume that 1 — ps a [s, h near ] -C 1 
for all S Q .] We can now state the new average-mismatch 
prescription as 



1 - ^(X'JAIAI 3 > MM, 



(87) 



which ensures that the detection rate, averaged over the 
extrinsic parameters of the signal, is reduced at most by 

3 

the factor MM . We shall call MM the average minimum 
match. 



B. Null parameter directions and reduced metric 

As discussed by Sathyaprakash and Schutz jiBJ and 
by Cutler 46], an extreme example of boundary effects 
occurs when one of the eigenvalues of <7bc a t X A (say, 
A(!)) becomes so small that it is possible to move away in 
parameter space along the corresponding eigendirection 
(say, efi ) and reach the boundary of the allowed param- 



eter region while keeping the mismatch S(X , X y 



i ) 



■(i) 



15 All the expressions to follow can be adapted to the case of a pri- 
ori known probability distribution for the extrinsic parameters. 
However, in our case it seems quite right to assume that the ori- 
entation angles © and <p are distributed uniformly over a sphere, 
and that a is distributed uniformly in the interval [0,27r]. 



well below the required value 5mm = 1 — MM. In other 
words, the ellipsoid of constant mismatch 5mm extends 
far beyond the allowed parameter region in the quasinull- 
eigenvalue direction. In such a situation, Eq. (|75|) will un- 
derestimate the total number of templates, because the 
denominator should now express the volume of the in- 
tersection of each lattice cell with the allowed parameter 
region. 16 A simple-minded fix to Eq. I|75|) is the follow- 
ing: write det gsc — Yi(k) ^(k)> where the Aqq are the 
n eigenvalues of gBC\ identify all the small eigenvalues, 
where small can be defined by A^) <C (1 — MM)/i^, 
with Z(j) the coordinate diameter of the allowed param- 
eter range along the eigenvector e A y, replace the small 
eigenvalues by the corresponding value of the expression 
(1-MM)//^; use this modified determinant in Eq. J7SJ). 

Physically, the presence of k small eigenvectors sug- 
gests that the variety of waveform shapes spanned by 
an n-dimensional template family can be approximated 
with very high overlap by an (n— fc)-dim. reduced family. 
A lower-dimensional template bank is certainly desirable 
for practical purposes, but it is necessary to exercise cau- 
tion: because the metric <7bc is not homogeneous, the 
quasinull eigendirections rotate as we move in parameter 
space, 1 so we need to show explicitly that any signal in 
the n-dim. family can be reached from a given (n — k)- 
dim. submanifold along a quasinull trajectory. For this to 
happen, the small eigenvalues must exist throughout the 
entire n-dim. parameter space, and the flow of the quasin- 
ull eigenvectors must map the submanifold into the entire 
space. To see that under these conditions the mismatch 
between the points on the submanifold and the points 
outside is indeed small, consider the following argument, 
due to Curt Cutler |4g. The triangle inequality for the 
inner-product distance guarantees that 



S 1/2 [X A (0),X A {1)} < 



dX B dX c 



9bc- 



dv dv 



-dv 



(88) 



along any path X A (v); for a path that follows the flow 
of the quasinull eigenvector eh (a reduction curve), the 
total mismatch is then bounded by the average of A/^ 
along the curve, times an integrated squared parameter 
length of order l(i)- 18 



16 Pictorially, the error that we make with Eq. IY5I is to let the 
template bank be thinner than a single template in the direction 

e A 

In fact, in the context of our templates this rotation is such that 
Eq. 1731 ceases to be true in the quasinull eigendirections for 
5 > 0.01. As soon as we move away from the point X A where 
the metric is computed, any rotation of the eigenvectors means 
that the original quasinull direction is no longer the path along 
which the mismatch grows most slowly. If the larger eigenvalues 
are several orders of magnitude larger than the smaller ones, 
as is true in our case, a tiny rotation is enough to mask the 
contribution from the smallest eigenvalue. 
18 At least if the geometry of the reduction curve is not very con- 
voluted. 
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Xi 

FIG. 8: Plot of (xi, ki) reduction curves in the (xi , fti) plane. 
We show curves for two sets of starting extrinsic parame- 
ters, corresponding to detector directions perpendicular (dark 
dots) and parallel (light dots) to the initial orbital plane. The 
curves start at the points marked with circles, and proceed 
in steps of 10 -6 for the nominal mismatch (i. e., the mis- 
match computed using the projected metric). For starting 
points at xi — 0.5, we follow the quasinull eigenvector for 
both positive and negative increments. The curves end at 
the (xi> K i) boundary, or (roughly) where the true mismatch 
(i. e., the exact mismatch between the local and the starting 
template) becomes greater than 0.01. The ending points are 
marked with crosses, and they are annotated with the num- 
ber of steps taken since the starting point, and with the true 
mismatch in units of 10 -3 . 



• 1st set of extrinsic parameters 




M(M e ) 

FIG. 9: Plot of (xi, ki) reduction curves in the (M, M) plane. 
The curves are the same as shown in Fig. |HJ but we omit all 
markings. 



For the ST^r template bank and for the two-stage 
search scheme of Sec.0 we find that the projected met- 
ric gfj° 3 admits a small eigenvalue for all values of the 
intrinsic and extrinsic parameters. Figures [S] and [5] show 
several examples of reduction curves that follow the quas- 
inull eigendirections [the subtleties related to projected- 
metric reduction curves are discussed in App. [Bj. The 
curves shown 19 begin at the points marked with circles, 
where (mi + m-i) = (10 + 1.4) M@ and 

(Xi,«i) = {5;jj}x | 0.0 J; (89) 

the curves then proceed in steps of 10 -6 for the nominal 
mismatch (i. e., the mismatch computed using the local 
projected metric) until they reach the (xi,ki) bound- 
ary, or (roughly) until the true mismatch (i. e., the exact 
mismatch between the local and the starting template) 
is greater than 0.01. We show curves for two sets of 
starting extrinsic parameters, corresponding to detector 
directions perpendicular (dark dots) and parallel (light 
dots) to the initial orbital plane. Figure [S] shows the 
projection of the reduction curves in the (xi, K\) plane; 
the ending points are marked with crosses, and they are 
annotated with the number of steps taken since the start- 
ing point, and with the true mismatch in units of 10~ 3 . 
Comparing the two numbers at each cross, we see that 
the triangle inequality is always respected: the true mis- 
match Sn is always less than the accumulated nominal 
mismatch 10 _6 iV 2 (where N is the number of steps); in 
fact, we see that the latter is a good approximation for 
the former. Figure El shows the projection of the same 
reduction curves in the (M,Ai) plane. The chirp mass 
M. = A/77 3 / 5 varies by less than 2% along the curves: this 
is natural, since M. dominates the evolution of the GW 
phase [see Eq. (0]. 

Figure |H1 suggests that we can reduce the dimensional- 
ity of our template bank by collapsing the (xi, K4) plane 
into ~ three curves, while retaining the full (M, rj) plane. 
Templates laid down on these 3-D submanifolds with a 
required minimum match MM will then cover every sig- 
nal in the full 4-D family with mismatch no larger than 



The curves of Figs. [5] and [5] are in fact obtained by following 
the quasinull eigenvectors of the fully projected (xi,^l) met- 
ric, which is a 2-D metric on xi and «i obtained by projecting 
9ij° 3 again over M and r\ [using Eq. 1831 1. as if they were ex- 
trinsic parameters. The two projection steps are equivalent to 
projecting the 9-D full metric into the 2-D (xi,^i) plane in a 
single step. This procedure estimates correctly the reduction 
mismatch introduced by adopting the reduced template family 
created by first collapsing the 2-D (xi,m) plane into several 1- 
D curves, and then including all the values of (M, n) for each 
point on the curves. We choose to collapse the (xi, Ki) plane for 
empirical reasons: the xi and k± parameter bounds are simple, 
(Xii K i) £ [Oil] x [ — 1 , 1] , and the reduction curves have large 
parameter lengths. 
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(1 — MM) + (5 rc( j, where <5 re d — 0.01 is the reduction mis- 
match introduced by the reduction procedure. Further 
investigations will be needed to find the optimal choice of 
reduction curves in the (xi, ki) plane, and to investigate 

the reduction curves of the average metric <?f 7 10J . 



C. Template counting 

While three or more reduction curves will probably be 
necessary to limit 5 ro( j — 0.01, for the sake of dcfinitcncss 
we select a 3-D reduced template space corresponding 
to (mi,m 2 ) S [1,3] x [7,12], ki = 0, and xi € (0, l]. 20 
We compute the total number of templates in this 3-D 
template bank according to 



templates 



jJ\detgp\dMdvd X 



(90) 



2y/(l - MM)/3 

where the primed indices i', j' run through M, rj, and Xi 

and we set A 4 = k\ — 0; furthermore, 3^ r ° J denotes the 
metric averaged over the extrinsic parameters 0, (p, and 
a, as given by Eq. (|86|l . The integral is carried out by 
evaluating the projected metric at the parameter sets 



(mi,m 2 ,xt) 



7M 
12M 







1M ' 






2M Q 


H 




3M Q 




0.7 




, 1.0 , 



at each of the points the metric is averaged on 100 pseu- 
dorandom sets of extrinsic parameters. The integra- 
tion then proceeds by interpolating across the param- 
eter sets (|9"TJ) . The final result is A/templates — 76,000 
for MM = 0.98 (not including the reduction mismatch). 
Given the uncertainties implicit in the numerical compu- 
tation of the metric, in the interpolation, in the choice 
of the reduction curves, and in the actual placement of 
the templates in the bank, this number should be un- 
derstood as an order-of-magnitude estimate. Most of the 
templates, by a factor of about ten to one, come from the 
parameter region near 777,2 = 1 (that is, from the small-ry 
region) . 



VII. SUMMARY 

Buonanno, Chen, and Vallisneri recently proposed 
[BCV2] a family of physical templates that can be used 



20 In fact, a second small eigenvector appears as we get close to 
Xi = 0; this is because spin effects vanish in that limit, so a 2-D 
family of nonspinning waveforms should be sufficient to fit all 
signals with small xi ■ 



to detect the GWs emitted by single-spin precessing bi- 
naries. The attribute physical refers to the fact that the 
templates are exact within the approximations used to 
write the PN equations that rule the adiabatic evolution 
of the binary. In this paper, after reviewing the definition 
of this template family (here denoted as ST^v), we discuss 
the range of binary masses for which the templates can 
be considered accurate, and examine the effects of higher- 
order PN corrections, including quadrupole-monopole in- 
teractions. We then describe an optimized two-stage de- 
tection scheme that employs the STjy family, and inves- 
tigate its false-alarm statistics. Finally, we estimate the 
number of templates needed in a GW search with LIGO- 
I. Our results can be summarized as follows. 

We determine the range of binary masses where the 
STn templates can be considered accurate by impos- 
ing two conditions: first, for the orbital separations 
that correspond to GWs in the frequency band of good 
interferometer sensitivity, the dynamics of the binary 
must be described faithfully by an adiabatic sequence 
of quasi-spherical orbits; second, the nonspinning body 
must be light enough that its spin will be negligible for 
purely dimensional reasons. The selected mass range is 
(mi, ma) ^ [7,12]M x [1,3]M . 

To evaluate the effect of higher-order PN corrections 
for binaries in this mass range, we compute the overlaps 
between templates computed at successive PN orders. 
When computed between templates with the same pa- 
rameters, such overlaps can be rather low; however, they 
become very high when maximized over the parameters 
(both intrinsic and extrinsic) of the lower-order PN tem- 
plate [see Table |n] . This means that the ST2 template 
family should be considered acceptable for the purpose 
of GW detection; but this means also that the estimation 
of certain combinations of binary parameters can be af- 
fected by large systematic errors [l9j • [When precessing- 
binary gravitational waveforms computed within PN- 
resummed and nonadiabatic approaches become 
available, it will be interesting to compare them with 
the PN-expanded, adiabatic STjv templates, to see if the 
maximized overlaps remain high. We do expect this to 
be the case, because the spin and directional parameters 
of the STtv templates provide much leeway to compen- 
sate for nontrivial variations in the PN phasing.] Again 
by considering maximized overlaps, we establish that 
quadrupole-monopole effects |40l |41| can be safely ne- 
glected for the range of masses investigated [Table ITTTj . 

We describe a two-stage GW detection scheme that 
employs a discrete bank of ST2 templates laid down 
along the intrinsic parameters (M, rj, Xij ^i) [although 
the (xij^i) ma y be collapsed to one or few 1-D curves, 
in light of the discussion of dimensional reduction of Sec. 
IVIj . The detection statistic ps<* (M, i],Xi> K i) is t ne over- 
lap between the template and the detector output, maxi- 
mized over template extrinsic parameters: (to, $0j Pi) = 
(tQ,$>o,9,(j),ip,@,<p). This maximization is performed 
semialgebraically, in two stages. First, for all possi- 
ble times of arrival to, we maximize the overlap over 
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$o and over Pj without accounting for the constraints 
that express the functional dependence of the Pj on 
(6, (f>, -0, ©, </?) : this step yields the approximated (uncon- 
strained) maximum p' Ea , which can be computed very 
rapidly, and which sets an upper bound for p-a . Sec- 
ond, only for the times of arrival to at which p' Ba passes 
the detection threshold, we compute the fully constrained 
maximum p-=a , which is more expensive to compute. 
[Note that this scheme differs from traditional hierarchi- 
cal schemes because we use the same threshold in the 
first and second stages.] We find that p' 3a is a good ap- 
proximation to /Ctea , so the number of first-stage triggers 
passed to the second stage is small. 

For a total false-alarm probability of 10 _3 /year, and 
for a conservative estimate for the number of independent 
statistical tests, the detection threshold is around 10. For 
this value, between 5 and 15 first-stage triggers are passed 
to the second stage for each eventual detection. For the 
same threshold, the single-test false-alarm probability is 
lower for ST2 templates than for the (^0 ^3/2^)6 DTF 
of BCV2 [the total false-alarm probability depends on 
the number of independent statistical tests, which is not 
available at this time for the (^0^3/2^)6 DTF]. 

The procedure of maximization over the extrinsic pa- 
rameters outlined in this paper can also be adapted for 
the task of detecting GWs from extreme-mass-ratio in- 
spirals (i. e., the inspiral of solar-mass compact objects 
into the supermassive BHs at the center of galaxies [471 ]') 
and inspirals of two supermassive black holes with LISA 
|48j . This is possible under the simplifying assumptions 
of coherent matched filtering over times short enough 
that the LISA antenna patterns can be considered con- 
stant, and of GW emission described by the quadrupole 
formula. Furthermore, the formalism of projected and 
reduced mismatch metrics developed in Sec. lVII can treat 
GW sources, such as extreme-mass-ratio inspirals, where 
many physical parameters are present, but only few of 
their combinations have significant effects on the emitted 
waveforms 45, 46]. In fact, this formalism is closely re- 
lated to the procedures and approximations used in the 
ongoing effort (motivated by mission-design considera- 
tions) to count the templates needed to detect extreme- 
mass-ratio inspirals with LISA |49| . 

It should be possible to generalize the formalism be- 
yond quadrupole GW emission, at least to some extent. 
When higher-multipole contributions are included, the 
detector response becomes much more complicated than 
Eq. (see, e. g., Eqs. (3.22b)-(3.22h) of Ref. [13). 
In particular, the response cannot be factorized into a 
factor that depends only on the dynamical evolution of 
the binary, and a factor that depends only on the po- 
sition and orientation of the detector; it is is instead a 
sum over a number of such terms, each containing dif- 
ferent harmonics of the orbital and modulation frequen- 
cies. Despite these complications, it should still be pos- 
sible to maximize the overlap over the extrinsic parame- 
ters, using a relatively small number of signal template 
and template-template inner products. The constrained- 



maximization procedure would however be very compli- 
cated, and although the (fully algebraic) unconstrained 
maximum would still be easy to compute, the dimension- 
ality of the unconstrained template space would now be 
so large that it may increase the false alarm probability 
too dramatically to make the two-stage scheme useful. 

The last result of this paper is an estimate of the num- 
ber of ST2 templates needed for a GW search in the mass 
range [7, 12] M© x [1,3]M Q . To obtain this estimate, we 
first compute the full mismatch metric, which describes 
the mismatch for small displacements in the intrinsic and 
extrinsic parameters; we then obtain the projected met- 
ric, which reproduces the effect of maximizing the overlap 
over the extrinsic parameters. At this point we observe 
that the projected metric has an eigenvector correspond- 
ing to a very small eigenvalue; this indicates that we 
can choose one of the four intrinsic parameters to be a 
function of the other three, so the dimensionality of the 
ST2 template bank can be reduced to three. For sim- 
plicity, we perform this reduction by setting ki = 0. We 
then compute the reduced mismatch metric, and obtain 
a rough estimate of ~ 76, 000 as the number of templates 
required for an average MM of 0.98, or 0.97 including an 
estimated reduction mismatch of 0.01. 
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APPENDIX A: ALGEBRAIC MAXIMIZATION 
OF THE OVERLAP OVER THE P 7 

In this section, we explore the algebraic maximization 
of p$ [see Eq. I|55fl]. given by 



P$ 



l A IJ P I P J 
B IJ PiPj 



(Al) 



over the Pj. We recall that the five Pj are combinations 
of trigonometric functions of three angles, and therefore 
must satisfy two constraints: luckily, both of these can 
be formulated algebraically. In light of the discussion of 
Sec. IIVBI the overall normalization of the Pi does not 
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affect the value of the overlap l|55fl. As a consequence, 
we can rescale the Pj and replace the first constraint by 

B IJ P I Pj = 1 , (A2) 

which enforces the normalization of the templates. This 
constraint is chosen only for convenience: the maximum, 
subject to this constraint, is exactly the same as the un- 
constrained maximum found by searching over the entire 
five-dimensional space. Let us work out its value, which 
will be useful later. Introducing the first Lagrangian mul- 
tiplier A, we impose 

-^ r [A IJ P I Pj-\(B IJ P I Pj-l)} = (A IJ -\B IJ )Pj = 0, 

(A3) 

which has solutions only for A corresponding to the eigen- 
values of AB _1 . For those solutions, we multiply Eq. 



(|A"3|) by Pjr to obtain 

X = A ij PjPj; (A4) 

using Eqs. (|A1|) and ljA2|) . we then see that A is the square 
of the overlap, so it should be chosen as the largest eigen- 
value of AB _1 . We then write the unconstrained maxi- 
mum as 



/c4° — max y/ max eigv A B 1 . (A5) 

By construction, p' Sa will always be larger than or equal 
to the constrained maximum, p-^a . 

The second constraint comes from Eq. (|58p. Writing 
out the STF components, we get 



'Pi 



det Pij = det —= 
v2 



-P 5 /\/3 

P2 

Pi 



-Pi 



P2 
+ P 5 /V3 
Pi 



P; 
Pa 

-2P 5 /V3) 



= D ijk PjPjP k = 0. 



(A6) 



[The tensor D IJK can be chosen to be symmetric since 
D UK P I P J P K = D^J^P^jPk] The constrained 
maximum of over the Pi, subject to the two con- 
straints, can be obtained as the maximum of the function 



A'-'PtPj - \(B 1J P!Pj - 1) - p(D 1JK PiP,,P K ) (A7) 

over Pi and over the two Lagrange multipliers A and /z. 
After taking partial derivatives, we get a system of seven 
equations, 



A IJ P, - XB 



"P., - °-pD IJK PjP K 



= 



B ij PiPj 
-.UK 



1 



D lJK PiPjP K = 



(A8) 

0, (A9) 
0, (A10) 



where the last two equations come from the con- 
straints (|A2I) and l|A6() . Multiplying the first equation 
by Pi and using the two constraints, we obtain Eq. (|A4() 
again. So the first Lagrange multiplier A is still the square 
of the overlap. The second Lagrange multiplier /i is zero 
when the signal s belongs to STn template family, and 
has exactly the same intrinsic parameters as the tem- 
plate. In this case, the extrinsic parameters of the signal 
correspond to a vector Pi that satisfies Eq. I|A8(I with 
p = (the multiplier A is still needed to enforce normal- 
ization of the template) . When the intrinsic parameters 
are not exactly equal, but close, (i becomes finite, but 
small. Equations (|A8|) - (jA10(l can then be solved itera- 



tively by expanding Pi in terms of fi, 



Pr 



(All) 



Inserting this expansion into Eqs. (|A8(I and (|A10|I . we get 

the zeroth-order equation 



/Jd(o) 



A u p 



(A LM Pi 0) P$)B IJ Pj 



(0) 



0. 



(A12) 



where we have already used the zeroth-order version of 
Eq. (|A4|> to eliminate A. 

Multiplying by (S _1 ) A/ , we see that the zeroth- 
order solution P^°' must lie along an eigenvector of 
(B~ 1 ) KI A IJ , and that the corresponding eigenvalue 

must be equal to A LM P^ pff , and therefore also to the 
square of the zeroth-order extremized overlap. To get the 
maximum overlap, we must therefore choose Pj ^ as the 
eigenvector corresponding to the largest eigenvalue. So 
the zeroth-order constrained maximum is exactly the un- 
constrained maximum obtained above [Eqs. (|A3|) - (|A5I) ]. 
We can then proceed to nth-order equations: 



[A JJ -2(A JM P$B IL P { ° } ) - (A LM P^'P^)B 1J }Py 

= E E ALMpi L 

mi=0 ni2— 
n-1 „ 



iLM p(0) p(0)x R /Ji p(») 



\LM p(mi) p("i 2 ) j^I J p(n-mi-m 2 ) _ 



E£ JjUK pM p(n-m-l) 
o J K 



(A13) 
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At each order, we insert the nth-order expansion of Pj 
into Eq. I|A10[) . and select the real solution closest to 
zero as the nth-order approximation to fi (such a solution 
is guaranteed to exist for all odd n). We then obtain 
the nth-order approximation to A (and therefore to ps<* ) 
using Eq. I|A4|I . We proceed in this way, until A and fi 
converge to our satisfaction. 

This iterative procedure succeeds when the intrinsic 
parameters of signal and template are close; as their dis- 
tance increases, the procedure becomes more and more 
unstable, and eventually fails to converge. The iteration 
fails often also when the overlap is optimized against pure 
noise. For these reasons, a practical implementation of 
the detection statistic ps° must eventually rely on the 
semialgebraic maximization procedure discussed in Sec. 
IIVBI Indeed, we have used the semialgebraic procedure 
for all the tests discussed in Sec. Ivl 



APPENDIX B: DIMENSIONAL REDUCTION 
WITH A NONUNIFORM PROJECTED METRIC 



In this appendix we extend the reasoning of Sec. IVI Bl 
to study dimensional reduction under the projected met- 
ric gfJ° 3 (^ A ), which lives in the intrinsic parameter space, 
but is a function of all parameters. For each point 
X A = (X' l ,E a ) in parameter space, we denote A(!)(A A ) 
the smallest eigenvalue of gfj° 3 (X A ), and e 1 ^ (\ A ) the cor- 
responding eigenvector in the intrinsic parameter space. 
Suppose we have 



A(i)(A' 4 )« 



1 - MM 
L (i) 



(Bl) 



for all values of X A in the allowed parameter region, where 
In) is the coordinate diameter of the allowed parameter 
range along the eigenvector . 

Now let us start from a generic point X A = (Xq, Sq ) in 
parameter space and follow the normal eigenvector 

for a tiny parameter length e, reaching X A = (X\,Ef), 
according to 

X[ = Xi + eefaiXf), (B2) 
SS + efr-W)] * 3 [C{X^)] p .e? (1) ^);(BS) 



— i 



this choice of AS Q makes the extrinsic parameter 
that minimizes 6(Xq, Sq ; X{, Sf). Denoting the inner- 
product distance as dist(A^, X A ) = a/2<J(Aq-, Xf), we can 
write 



dist(A^, X A ) = e^2A (1) (X A ) + 0(e 



(B4) 



from X A , we follow the eigenvector e^(A^) for another 

parameter length e, and reach X A ; then from X A we reach 
X A , and so on. Up to the Ath step, we have traveled a 




ductionP^ 



FIG. 10: Illustration of dimensional reduction. Here we show 
a signal space with one extrinsic parameter (H 1 ) and two in- 
trinsic parameters (X 1,2 ), and we assume that the projected 
metric has one small eigenvalue all through parameter space. 
Starting from a generic point Xq, we follow the flow of the 
quasinull eigenvector of gfj 01 for an infinitesimal parameter 
distance to reach Aj 4 ; we then repeat this process, each time 
adjusting the direction of the eigenvector according to the 
metric (hence the difference between the reduction path pre- 
dicted at Xq and the actual reduction path). In the end we 
reach Ajv after having accumulated a parameter length I in 
the intrinsic parameter space. The mismatch between Ao and 
Ajv will be smaller than 5mm = 1 — MM, if I is not much larger 
than In), the coordinate diameter of the intrinsic parameter 
space in the approximate direction of the quasinull eigenvec- 
tor. 



cumulative parameter distance I = Ne in the intrinsic 
parameter space, and an inner-product distance 



dist(A^,A^) < ^dist(Al 1 ,A^) 

n=l 

N r 

= e^2A^U)+0(e 2 ) 



< l^2m&xA {1) (X A ) + 0(Ne 2 ), (B5) 

where in the first line we have used the triangle inequal- 
ity for the inner-product distance. The term 0(Ne 2 ) 
vanishes in the limit e — * 0, N — > oo, keeping I = Ne 
finite (see Fig. I1U|) . So we can take the continuous limit 
of Eqs. (|B2(1 and (|B3(1 and arrive at two differential equa- 
tions that define the resulting trajectory: 



X i (l) = e\ 1) , 5 a (Z)=[ 7 - 1 ]^C w e| 1) 



(B6) 



where X 1 and S a arc parametrized by the cumulative 
parameter length I, with 



X\l = Q) = X l a , S Q (/ = 0) = S« 



(B7) 



We can allow I to be either positive or negative, in order 
to describe the two trajectories that initially propagate 
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full signal space x 




Y / & 




space 


X l 





FIG. 11: Illustration of reduced signal space as a hypersur- 
face inside full signal space. Here we show only the direc- 
tions along the intrinsic parameters. Starting from the points 
(X[,3"), (XI, S"), we follow the trajectory HB6H and 

reach the hypersurface at (X{ red , S? red ), . . . , (X\ xed , E£ red ). 
For these particular initial points, Xi re< j happens to coincide 
with Xirod, and X2red with Xg Ici i- We can see that and 
\£ (and indeed all points that reduce to X\ rc d, including the 
points along the reduction curve) will be indistinguishable 
upon detection with the reduced template bank. The same is 
true for A^, A3 1 , and for all the points that reduce to -Xi le d- 



along ±e^(Ao )■ Equation (|Bf>|) then becomes 



dist [A^,A A (0] 



< 



< 



dl'^2A {1) [\ A (V)\ 
|/|^2maxA (1) (A A ). (B8) 



In terms of mismatch, 



min<5 [X^-X^l)] = ^ mm dist [\ A ;X i (l)] 1 



< 



< 



lf dist [A^; A a (/)] 



dl' y/2A (1) [X A (l')} 



< Z 2 maxA (1) (A A ) 



(i) 



~r~ I S 



J MM . 



(B9) 



where the hybrid notation of the first line indicates the 
mismatch along the solution of (|B6(I . and where of course 
8mm = 1 — MM. Here, although we evolve X 1 and E a 
simultaneously, it is the trajectory X l (l) in the intrinsic 
parameter space that we are ultimately interested in. In 
the context of dimensional reduction for the projected 
metric, we shall call A" 4 (7) the reduction curve. 

If the reduction curves are reasonably straight, it 
should be easy to find a (dimensionally reduced) hyper- 
surface with the property that any given point (Xq, Eq) 
in the full parameter space admits a reduction curve 
that reaches the hypersurface at a parameter 2* not 
much larger than the coordinate diameter of parame- 
ter space (see Fig. 1110 . From Eq. (|B9|) . we then have 
mins" S[Xq, Eq ; X l (l r )] < 5mm- So any point in the full 
parameter space can be fit with a mismatch smaller than 
<5mm by a point on the hypersurface. 
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